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6.  Abstract 

A three-dimensional  mathematical  model  of  an  aircraft  seat,  occupant,  and  restraint 
system  has  been  developed  as  an  aid  to  the  development  of  crashworthy  seats  and 
restraint  systems  for  general  aviation  aircraft.  The  occupant  model  consists  of 
eleven  rigid  mass  segments  whose  dimensions  and  inertial  properties  have  been 
determined  from  studies  of  human  body  anthropometry  and  kinematics.  The  seat  model 
is  made  up  of  beam  and  membrane  elements  with  provision  for  simulating  plastic 
behavior  by  the  introduction  of  plastic  hinges  in  the  beams. 

A user-oriented  computer  program'called  Seat  Occupant  Model-Light  Aircraft  (SOM-LA) 
based  on  the  three-dimensional  model  has  been  developed  for  use  by  engineers 
concerned  with  design  and  analysis  of  general  aviation  seats  and  restraint  systems 
in  that  detailed  descriptions  of  both  are  used  as  input.  The  response  of  the  seat 
and  occupant,  restraint  system  loads,  and  various  injury  criteria  are  predicted  for 
any  given  set  of  crash  conditions. 
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I Q INTRODUCTION 


The  design  of  crashworthy  seats  and  restraint  systems  for 
aircraft  presents  a complex  engineering  problem,  the  solution  of 
which  can  be  greatly  aided  by  sufficiently  rigorous  analytical 
techniques.  The  crash  environment  can  vary  widely  from  one  acci- 
dent to  another  such  that  a great  number  of  conditions  must  be 
evaluated  to  establish  those  critical  to  occupant  survival.  For 
example,  the  restraint  system  must  limit  the  movement  of  the  occu- 
pant sufficiently  to  eliminate  the  possibility  of  head  strike  on 
rigid  cockpit  structure.  Also,  the  relatively  low  tolerance  of 
the  human  body  to  accelerations  in  a direction  parallel  to  the 
spine  requires  the  consideration  of  vertical  impact  forces 
which  are  usually  present  and  often  significant  in  aircraft 
crashes.  It  is  obvious  that  a very  strong  seat  is  not  a valid 
solution  for  this  problem,  since  it  would  not  only  incur  serious 
weight  penalties  but  would  transmit  high  vertical  impact  forces 
directly  to  the  occupant.  Rather,  the  design  of  a crashworthy 
seat  for  aircraft  includes  the  capacity  to  absorb  energy  through 
controlled  deformation  in  the  vertical  direction,  thus  lowering 
the  accompanying  loads. 

In  the  initial  phases  it  is  desirable  to  evaluate,  in  some 
detail,  existing  seats  and  restraint  systems  in  their  surrounding 
cockpits,  thus  establishing  existing  weaknesses.  It  is  then 
desirable  to  make  modifications  and  to  evaluate  the  effect  of 
these  modifications  on  improving  the  survivability  of  the  system. 
These  evaluations  must  be  conducted  for  a great  many  of  the  pos- 
sible crash  environments,  thus  constituting  a relatively  large 
matrix.  It  is  apparent  that  testing  would  be  extremely  expen- 
sive and  require  a great  deal  of  time,  since  design  modif icat:'  ,-ns 
would  have  to  be  developed  and  fabricated  prior  to  testing. 
Therefore,  an  analytical  technique,  such  as  was  developed  in 
this  program,  is  required, 

A number  of  one-,  two-,  and  three-dimensional  mathematical 
models  of  the  human  body  have  been  developed  for  crash  surviva- 
bility analysis.  These  simulation  models  have  generally  been 
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intended  for  use  in  evaluation  of  automobile  interior  design  with 
respect  to  injuries  caused  by  secondary  impacts,  such  as  the 
three-dimensional  models  described  in  references  1 through  3. 

Seats  have  been  represented  in  a very  simple  manner  because  in 
automobiles  the  role  of  the  seat  design  in  determining  occupant 
survival  is  minimal.  A simulation  model  intended  specifically 
for  this  application  was,  therefore,  required. 

The  development  of  a three-dimensional  mathematical  model  of 
a light  aircraft  seat,  occupant,  and  restraint  system  is  described 
in  this  report.  This  model  forms  the  basis  for  a simulation  com- 
puter program  that  has  been  written  specifically  for  use  in  crash- 
worthy  design  and  analysis  of  light  aircraft  seats  and  restraint 
systems.  The  program  has  been  organized  so  as  to  minimize  the 
volume  and  complexity  of  input  data  and  to  focus  on  seat  and  re- 
straint system  design  parameters.  The  effort  described  herein 
was  performed  for  the  Federal  Aviation  Administration,  Systems 
Research  and  Development  Service,  under  Contract  DOT-FA72WA-3101, 
which  commenced  in  August  1972. 
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2.0  MATHEMATICAL  MODEL 

The  three-dimensional  mathematical  model  includes  a lumped- 
parameter  representation  of  the  vehicle  occupant  and  a finite  ele- 
ment seat.  Interface  between  the  seat  and  occupant  is  provided 
by  seat  cushions  and  a restraint  system,  consisting  of  a lap  belt 
and,  if  desired,  a single-strap  or  double-strap  shoulder  harness. 
The  response  of  the  occupant  and  seat  can  be  predicted  for  any 
given  set  of  aircraft  impact  conditions,  including  the  initial 
velocity  and  attitude  and  the  input  acceleration. 

This  section  provides  a discussion  of  the  development  of  the 
occupant  and  seat  models,  including  details  of  the  approach  to 
formulating  the  equations  of  motion  and  or  the  technique  used  for 
their  solution.  The  first  part  of  this  section  discusses  the  de- 
velopment of  the  equations,  and  subsequent  parts  cover  particular 
aspects  of  the  equations,  specifically  the  body  joint  model  and 
the  treatment  of  external  forces. 

2.1  OCCUPANT  MODEL 

The  mathematical  model  of  the  aircraft  occupant  is  made  up 
of  11  rigid  segments,  as  shown  in  figure  1.  This  number  is 
thought  to  represent  the  minimum  that  will  permit  accurate,  mean- 
ingful simulation.  A greater  number  might  possibly  improve  the 
accuracy  of  simulation  but  would,  in  turn,  increase  program  exe- 
cution cost.  Arm  and  leg  segments  are  included  to  enable  predic- 
tion of  injuries  to  these  extremities.  Although  leg  and  arm 
injuries,  in  themselves,  may  not  be  as  serious  as  head  or  chest 
injuries,  they  may  prevent  escape  from  a stricken  aircraft  and  the 
potential  hazard  of  post.crash  fire. 

Each  of  the  body  joints,  with  the  exception  of  the  elbow  and 
knee  joints,  possesses  three  rotational  degrees  of  freedom.  Be- 
cause of  the  hinge-type  motion  of  a forearm  or  lower  leg  relative 
ro  an  upper  arm  or  thigh,  respectively,  the  position  of  each  of 
these  segment-:  is  described  by  one  additional  angular  coordinate. 
Therefore,  the  occupant  system  possesses  a total  of  28  degrees  of 
f) aedom. 
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2.1.1  Equations  of  Motion 

2. 1.1.1  Coordinate  Systems 

Fixed  at  the  center  of  mass  of  each  of  the  11  segments  is 
a right-handed  cartesian  coordinate  system.  For  segment  n (n  = 1, 
2,  . ..,  11)  the  local  coordinate  system  is  denoted  by  axes  (x  , 
y , 2n) . Positive  directions  are  defined  such  that  when  the  body 
is  seated  as  shown  in  figure  2,  with  the  torso  and  head  upright, 
the  upper  arms  parallel  to  the  torso,  and  the  elbows  and  knees 
bent  at  right  angles,  xn  is  positive  forward,  y^  to  the  left,  and 
zn  upward. 
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In  order  to  describe  a general  position  of  the  body,  it  is 
necessary  to  relate  the  orientation  of  each  segment  (xR,  y^,  zR) 
to  the  inertial  system  (X,  Y,  Z) . The  angular  relationship  be- 
tween the  local,  segment-fixed  coordinate  and  the  inertial  system 
can  be  expressed  by  the  transformation 


(1) 


Because  three  angular  coordinates  can  be  used  to  define  the 
rotation  of  a given  segment,  it  is  convenient  to  utilize  a set  of 
coordinates  that  will  suffice  as  generalized  coordinates  in  the 
formulation  of  the  equations  of  motion.  A system  of  Euler ian 
angles  provides  a convenient  set  of  three  independent  angular 
coordinates.  Assuming  that  the  local  (xr,  yR,  z ) system  is  ini- 
tially coincident  with  the  inertial  (X,  Y,  Z)  system,  the  Euler 
angles  are  a series  of  three  rotations,  which,  when  performed  in 
the  proper  sequence,  permit  the  system  to  attain  any  orientation 
and  uniquely  define  that  position.  The  particular  set  of  Euler 
angles  selected  for  use  here  is  illustrated  in  figure  3 and  de- 
fined, as  follows: 

1.  A positive  rotation  \p  about  the  Z-axis,  resulting  in 
the  primed  (x',  y',  z'}  system. 

2.  A positive  rotation  9 about  the  y'-axis  resulting  in 
the  double-primed  (x",  y",  z")  system. 

3.  A positive  rotation  <p  about  the  x"-axis  resulting  in 
the  final  (x,  y,  z)  system. 

In  order  to  determine  the  elements  of  the  transformation 
matrix  [Tn] , it  is  necessary  to  consider  the  matrix  equations 
that  indicate  the  three  individual  rotations  described  above. 
Referring  again  to  these  definitions  of  0,  and  <j>,  the  following 
equations  are  obtained: 
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Writing  equations  (2)  - (4)  in  abbreviated  form 


{R}  = M {r*} 

{r ' } « [0]  {r"} 
(r" } = I4>3  {r} 

or 

(R>  = W [63  14>]  {r} 


where  {R}  represents  the  components  of  a vector  in  the  inertial 
system  and  {r}  represents  the  same  vector  in  the  final  (x,  y,  z) 
system.  Performing  the  matrix  multiplications  indicated  in 
equation  (5) , the  elements  of  the  transformation  matrix  in 
equation  (1)  are  obtained: 
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T..11  = COS  ll)  cos  0 

11  n n 

T.  _n  = cos  ij > sin  9 sin  <j>  - sin  ill  cos  4> 

12  n n n n n 

T.  * cos  tj>  sin  0„  cos  <b  + sin  ijj  sin  <|> 

13  n n n n n 

T0 1 n = sin  to  cos 

21  n n 

T00n  = sin  to  sin  0„  sin  to  + cos  to  cos  to 

22  Tn  n n n n 

Tnr>n  = sin  to  sin  0 cos  to  - cos  to  sin  to 

..3  n n n n n 

T,.n  = -sin  0 
3j.  n 

T,,.n  = cos  0 sin  <J> 

32  n n 
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(6 


= cos  0_  cos  4 


for  n = 1,  2,  3,  4,  6,  8,  10 

The  additional  constraint  of  hinge-type  rotation,  at  the 
elbows  and  knees,  requires  the  use  of  one  additional  angular 
coordinate  to  define  the  position  of  each  of  the  forearm  and 
lower  leg  segments.  Referring  to  figure  4,  the  angular  position 
of  the  forearm  segments  (1-5,  7)  is  given  by 


Figure  4.  Definition  of  Angular  Coordinates  a for  Elbows 
and  Knees. 
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and  the  lower  leg  segments  (m  = 9,  11)  by 


m 


m 


m 


sin  a 0 -cos  a 
m m 


0 10 
cos  a 0 sin  a, 


m 


m 


m 


‘m 


m 


(8) 


From  equations  (6)  and  (7)  the  elements  of  the  transformation 
matrix  for  the  forearms  (sements  5 and  7)  are  written  as: 

0 

T, , = cos  to  cos  6 sm  ct„  - cos  to  sin  6..  cos  <j>  cos  a„ 

11  n n £ n n n £ 


- sin  \p  sin  4>  cos  a. 
n n £ 


T1  - - cos  to  sin  0„  sin  to  - sin  to  cos  to 

12  n n n n n 


T13  = cos  ^n  cos  6n  cos  a£  + cos  ^n  sin  0n  sin  ^n  sin  aJ 


+ sin  sin  to  sin  a. 

n n £ 


T-,  = sin  to  cos  0 sin  a„  - sin  to  sin  0„  cos  to  cos  a0 

21  n n £ n n n £ 


- cos  to  sin  <j>  cos  a„ 
n n x, 


= sin  sin  0 sin  to  + cos  to  cos  4> 

22  n n n n n 

£ I 

= sin  to  cos  0 cos  a,  + sin  4>  i sin  0„  cos  to  sin  a, 

23  n n £ n n n J 


- cos  to  sin  to  sin  a0 
n n £ 


T 


31 


- sin  0 sin  a.  - cos  0 cos  to  cos  a„ 
n £ n n £ 


32 


cos  0 sin  4> 
n Yn 


'33 


- sin  0 cos  a„  + cos  0 cos  <j>  sin  a„ 
n £ n n £ 


(9) 
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From  equations  (6)  and  (8)  the  elements  of  the  transformation 
matrix  for  the  legs  (segments  9 and  11)  are  obtained: 


T. . = cos  cos  0 sin  a + cos  A sin  0_  cos  A cos  a 


+ sin  A sin  <p  cos  a 
n n m 


T.  = cos  A sin  sin  <p  - sin  A cos  A 

■ y 9 ii  9 n * w • !• 


T _m  = cos  ip  cos  6 cos  a + cos  A sin  0„  cos  A sin  a 


+ sin  ip  sin<p  „ sin  a 
n n m 


T01"  = sin  A cos  9^  sin  a + sin  A sin  6^  cos  A cos  a 

J I Tl  Ti  TH  ' Tl  n 1 m Ti 


- cos  A sin  A cos  a 
n n m 


T__  = sin  A_  sin  9 sin  A + cos  A cos  A 


T0  n = - sin  A„  cos  0,,  cos  a + sin  A sin  0„  cos  A sin  a 

/ A Tl  Tl  m n Tl  Tl  n 


- cos  A sin  ip  sin  a 
n n m 


sin  0„  sin  a + cos  0 cos  A cos  a 
n m n n m 


cos  0 sin  A 
n Yn 


sin  0„  cos  a + cos  0„  cos  <p  sin  a 
n m n n m 


Having  developed  the  relationships  expressed  in  equations 
(1)  through  (10) , the  position  of  the  occupant  can  be  described 
by  the  following  set  of  generalized  coordinates: 


*1 

*2 

q3 

q4 

q5 

<*6 

^7 

q8 

q9 

q10 


= X, 


* Z1 
= 

- 0i 


= <j> 


1 
*2 
02 
*2 
h 


qn 

q12 

q13 

q14 

q15 

q16 

q17 

q18 

q19 

q20 


* 63 
- *3 
= n 

= 04 

= *4 

= a5 

= *6 

= 96 

= $6 

= a-. 


q21 

q22 

q23 

q24 

q25 

q26 

q27 

q28 


= 4> 


8 


= 0 


8 

*8 


= a. 


= 'P 
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* 0 


10 

*10 

all 


(11) 


The  above  coordinates  include  the  cartesian  coordinates  of  the 
mass  center  of  segment  1 (X^,  Y^,  Z^) , selected  as  a reference 
point  on  the  body,  seven  sets  of  Eulerian  angles,  and  the  four 
additional  angular  coordinates  for  the  elbows  and  knees. 

2. 1.1. 2 Lagrange’s  Equations 

The  response  of  the  occupant  system  is  described  by 
Lagrange's  equations  of  motion,  which  are  written  for  the  28 
generalized  coordinates.  The  equations  are  developed  according 
to 


d ,3L  . 

dt  (7^') 

3qj 


= Q, 

3qj  3 


(j  = 1,  2,  ...,  28) 


(12) 


12 


where  L is  the  Lagrangian  function 


L = T - V 


(13) 


t represents  time.  Q.,  are  the  generalized  forces  not  derivable 
from  a potential  function » (Forces  that  are  derivable  from  a po- 
tential function  are  obtained  from  L and  T and  V are  the  system 
kinetic  and  potential  energies,-  respectively.) 

Because  the  system  being  treated  does  not  involve  any 
velocity-dependent  potentials,  equation  (12)  can  be  written  as 


d ,3T  . 3T  3V  n 

m aq,  ' + " 3 


(j  = 1,  2.  ....  28) 


(14) 


The  system  kinetic  energy  contains  both  translational  and  rota- 
tional parts: 


11 


T - 7 E M„  ux.)2  * (t;.?  T up2! 

n=l 


11 

+ j E (I 

n=l 


2 . 


x x 
n n 


y*  yn 


X z 
n n 


where  M is  the  mass  of  segment  n and  I - I , and  i as  mas*? 

n -v  Y‘  ” 

n - n n 

moments  of  inertia  of  segment  n with  respect  to  the  ]^ca?  coordi- 
nate aces  (x  , y , z ) , assumed  to  be  principal  moments  of  iner- 
n Jn  n ' r 

tia. 

The  absolute  velocities  of  the  11  mass  segments  required 
for  the  translational  kinetic  energy  must,  of  course,  be  written 
as  functions  of  the  generalized  coordinates  and  generalized  veloc- 
ities in  order  to  use  equation  (14) . The  derivation  of  these 
velocities  is  presented  in  appendix  A.  The  angular  velocity  com- 
ponents (w  , w , u)  ) seen  in  equation  (15)  are  parallel  to 
X V z 

n n 

the  local  (x  , y , zn)  coordinate  systems.  These  angular  veloc- 
ity components  cannot  be  used  directly  in  Lagrange’s  equations 
because  they  do  not  correspond  to  the  time  derivatives  of  any 
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set  of  coordinates  that  specify  the  position  of  the  segment. 

They  must  be  written  as  functions  of  the  generalized  coordinates, 
using  the  generalized  angular  velocites  (i^,  8^,  , which  are 

parallel  to  the  axes  Z,  y^,  and  x£,  respectively. 


An  arbitrary  angular  velocity  of  segment  n,  wn, 


can  be  ex- 
pressed as  a function  of  the  generalized  angular  velocities 
according  to 


co  = li)  + 

~n  -n 


+ <p 

-n 


(16) 


Referring  to  figure  3,^,8,  and  $ do  not,  in  general,  form  a 

— n — n — n 

mutually  perpendicular  vector  triad.  ($  and  | are  both  perpen- 
dicular to  Q but  are  not  necessarily  perpendicular  to  each  other.) 
However,  they  can  be  considered  as  a nonorthogonal  set  of  compo- 
nents of  w since  their  vector  sum  is  equal  to  u.  Summing  the 
orthogonal  projections  of  $n,  Q_n,  and  on  the  (xn,  yR,  yR) 
axes  yields  the  angular  velocity  components  required  for  the 
kinetic  energy  expression: 


- iL  sin  6 
Tn  n 

cos  sin  + 
n n 

cos  8 cos  <J) 
n n 


0n  cos  4,n 
en  sin  4,n 


(17) 


All  quantities  required  to  develop  the  system  kinetic  energy 
are  now  available,  and  the  entire  expression  is  written  in  appen- 
dix B. 


The  system  potential  energy  is  simply  gravitational  poten- 
tial, which  is  written  as 


11 

v=  £ 

n=l 


M g 


(Z  - 
n 


n 


(18) 
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where  g is  the  acceleration  due  to  gravity  and  ZR  is  an  arbitrary 

o 

datum.  The  potential  energy  expression  is  expanded  in  appendix  C. 
2 . 1 . 1 . 3 Matrix  Equations 

For  purposes  of  computation/  the  equations  of  motion  are  re- 
written in  the  following  form: 


[A(q)]  {q}  = (B(q,q)}  + {P  (q)  > + (R(q,q)}  + (Q(q,q)>  (19) 

where  the  elements  of  the  ine.rhia  matrix  [A]  and  the  vector  {b} 
are  derived  from  the  kinetic  energy  derivatives  of  Lagrange's 
equations.  In  other  words, 


d_ 

dt 


3T 


28 

E 

k=l 


AjA 


“ •••#  ^28'  ^1'  ^28^ 


(j  = 1,  2,  28)  (20) 

The  elements  of  [A]  and  (B } are  given  in  appendix  D.  The  force 
vector  { P } is  derived  from  the  system  potential  energy  according 
to 


P (^li*  • * * i q)  ~ ~ ^ 3 = if  2 , « « • , 2G) 

3 -L  aqj  (21) 

The  elements  of  { P } are  presented  in  appendix  E.  Both  {R}  and 
{ Q } are  vectors  of  generalized  forces  derived  from  the  right-hand 
side  of  Lagrange's  equations.  The  vector  (R } describes  the  re- 
sistance of  the  body  joints  to  rotation  and  is  discussed  in  detail 
in  section  2.1.2.  {Q}  is  the  vector  of  generalized  external 

forces  and  is  discussed  in  detail  in  section  2.1.3. 
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2.1.2  Joint  Resistance 

The  form  of  the  joint  resistance  vector  {R}  in  equation  (19) 
depends  on  the  user's  selection  of  occupant  - either  dummy  or 
human.  Although  both  joint  models  contain  the  same  types  of  ele- 
ments, a nonlinear  torsional  spring  and  a viscous  torsional 
damper,  the  function  of  each  of  these  elements  determines  the 
type  of  occupant. 

The  10  body  joints,  illustrated  in  figure  5,  are  defined  as 
follows : 

Joint  1 - Back,  between  12th  thoracic  and  1st  lumbar 

vertebrae 

Joint  2 - Neck,  between  7th  cervical  and  1st  thoracic 

vertebrae 

Joint  3 - Right  shoulder 

Joint  4 - Right  elbow 

Joint  5 - Left  shoulder 

Joint  6 - Left  elbow 

Joint  7 - Right  hip 

Joint  8 - Right  knee 

Joint  9 - Left  hip 

Joint  10  - Left  knee 

The  angular  displacement  of  joint  i from  its  reference  posi- 


( , / j 

m'  Jm 


fe  ) and  U 


m 


n' 


'n' 


k ) 

n 


tion  (figure  2)  is  given  by  B^.  If 
are  triads  of  unit  vectors  in  the  local  coordinate  systems  of  two 
adjacent  segments  connected  at  joint  i,  as  shown  in  figure  6,  the 
joint  angle  is  given  by 


B^  = cos 


-1 


(k 


m 


fe„> 


(22) 


where  • kn)  is  the  scalar  product.  Considering  the  geometry 
of  the  occupant  model  in  the  reference  position,  the  B^  for  the 
10  joints  are  given  by 
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m 


t 


Figure  6.  Joint  Angle  8^  Between  Segments  m and  n. 


B-l  = cos"1  (k1  * k2 ) 


B0  = cos"1  (k-  ’ fe,) 


B3  = cos"1  (fe2  * fe4) 


B4  - a5 


B5  = cos"1  (fe2  * feg) 


B6  = a7 
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! 


8?  = cos-1  * -tg) 

38  = “9 

gg  = cos  ~1  U x * -t1Q) 


3i0  = all 


(23) 


If  at  each  joint  i,  a momenL  and  a torsional  damper  with 
coefficient  act  to  resist  motion  of  the  joint,  then  the  virtual 
work  done  on  the  system  as  each  joint  i undergoes  a virtual  dis- 
placement 6 8^  is 


10 

6W  = - £ (M.  68.  + J.  B.  68.) 

i=l  11  ill  (24) 

Since  the  8^  are  functions  of  the  generalized  coordinates  q^, 
the  virtual  displacement!?.'  68^  can  be  expressed  in  terms  of  cor- 
responding virtual  displacements  of  the  q ^ . In  general,  such  an 
expression  would  take  the  form 


68±  = 


28  98,- 

jS  ^ Sq: 


(i=l , 2 , . . . , 10) 


(25) 


where  the  partial  derivatives  38^/3qj  are  functions  of  the  gen- 
eralized coordinates.  Substituting  into  equation  (24)  gives 


10 

SW  = - £ 


28 


38, 


L CM.  + J,  8j)  6q. 

i=l  j = l 1 11  9qj  3 


(26) 


Changing  the  order  of  summation,  equation  (26)  can  be  written  in 
the  general  form 


28 

£W  = £ Q. 

j-1  3 3 (27) 
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where  Q..  are  the  generalized  forces  acting  on  the  svstem.  As  seen 
in  equation  (19) , the  generalized  forces  are  being  treated  as  two 
distinct  types:  joint  resistance  forces  and  external  forces. 

Since  the  joint  resistance  terms  are  being  treated  here,  the  gen- 
eralized joint  forces  referred  to  as  R.  will  be  considered  alone. 
Equation  (27)  becomes,  more  specifically 


28 

6W  = E R-t  6q. 
j-1  3 3 


(28) 


and,  from  equation  (26) , can  be  written 


10 


4 = - E (M.  + J.  6.)  ffi  (j-1,  2,  ,28) 


i=l 


9q 


j 


(29) 


The  elements  of  vector  {R}  and  their  derivation,  using  equa- 
tion (29) , are  preserved  in  appendix  P. 

As  mentioned  earlier  in  this  section,  the  type  of  occupant 
is  determined  by  the  relative  contributions  of  and  to  the 
Rj  terms.  For  the  dummy  joint  the  resisting  torque  is  con- 
stant throughout  the  normal  range  of  joint  motion  and  increases 
rapidly  along  a third-order  curve  to  a higher  value  at  the  limit- 
ing displacement  , as  shown  in  figure  7.  The  normal  values 
Mp  are  set  equal  t&  those  resulting  from  the  joint-tightening 
procedure  of  SAE  Recommended  Practice,  Anthropomorphic  Test  De- 
vice for  Dynamic  Testing  - SAE  J963.  That  is,  the  body  joints 
will  just  support  a 1G  load  in  the  reference  (seated)  position, 
with  the  exception  of  the  torso  joints,  which  will  support  a 2G 
load.  In  addition  to  , a small  viscous  damping  term  with  con- 
stant is  included  for  energy  dissipation. 

The  resistance  of  the  human  joint  consists  of  up  to  three 
terms.  The  primary  resisting  force  during  normal  joint  rota- 
tion is  a viscous  damping  term  with  constant  coefficient  In 

a manner  similar  to  the  case  of  the  dummy,  a resisting  torque  is 
applied  at  the  limit  of  the  joint  range  of  motion,  as  shown  in 
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t 


ANGULAR  DISPLACEMENT  - $i 
Figure  7.  Dummy  Joint  Resisting  Torque. 


figure  8(a).  An  additional  term  used  to  simulate  muscle  tone  is 
the  moment  M* , which  drops  to  zero  after  a small  angular  dis- 
placement from  the  initial  position,  provided  that  the  crash  de- 
celeration is  sufficient  to  overcome  it  (figure  8(b)).  A com- 
pletely relaxed  occupant,  usually  the  worst  case  with  regard  to 
injury,  can  be  simulated  by  setting  this  last  joint  resistance 
term  to  zero. 


k 

f: 
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Figure  8.  Human  Joint  Resisting  Torques:  (a)  Displacement- 

Limiting  Moment;  (b)  Muscular  Resistance. 


2.1.3  External  Forces 

The  vector  of  generalized  external  forces  {Q}  is  developed 
in  a manner  similar  to  that  discussed  in  the  previous  section  for 
the  joint  resistance  vector.  The  resultant  external  force  r\ 
acting  on  segment  i is  given  by 


F.  = Fy  l + F / + F7  fe 
-x  X.  i 


(30) 


where  Fv  , Fv  , and  F„  are  components  in  the  inertial  (X,  Y,  Z) 
Xi  Yi  Zi 

system.  Let  the  absolute  position  of  the  point  on  segment  i, 
where  the  resultant  force  acts,  be  represented  by 


-P.  XP.  * + YP.  j + ZP. 

ii  i i (31) 


As  the  resultant  force  applied  to  each  segment  i undergoes  a vir- 
tual displacement  <$r  , having  components  (<$X0  , 6Y_  , 6Zn  ) , the 

i Pi  Pi  Pi 

virtual  work  on  the  system  done  by  the  F^  is 
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(32) 


SW=  £ (Px.  «p.  + F «Yp  + P St  ) 

1=1  11  11  11 

Writing  the  virtual  displacement  components  in  terms  of  the  gener- 
alized coordinates  q ^ : 


results  in 

28  11 

«W  = E E 

j=l  i=l 

Using  equation  (27) 
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9X 

p. 

6XP.  - 
1 

E 

j-i 

l 

3qj 

Sq. 

28 

3Yp. 

6YP.  = 
1 

E 

j-i 

i 

9qj 

28 

9V 

6Zp  = 

T_ 

E 

j=i 

i 

3qj 

fojj 

9X 

3YP 

(F 

xi 

Pi 

3qj 

+ F 

Yi 

(33) 


9Z 


+ F 


Z.  9q 
l ^3 


-)  <Sq^ 


(34) 


28 

«W  = I Q,  <5q^ 
j = l -1  J 


yields  the  components  of  the  generalized  external  force  vector: 


9X 


P. 

l 


9Y, 


9Z 


p. 


■i) 


(35) 


The  components  of  {Q}  derived  from  equation  (35)  are  presented  in 
appendix  G. 
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The  external  forces  acting  on  the  11  body  segments  can  be 
characterized  as  either  contact  forces  or  restraint  forces. 

These  forces  are  discussed  in  further  detail  in  the  sections 
following. 

2. 1.3.1  Contact  Forces 

The  contact  forces  exerted  on  the  occupant  include  all 
those,  except  the  restraint  forces.  As  illustrated  in  figure  9, 
these  forces  are  exerted  by  the  cushions,  floor,  and  an  optional 
air  bag  restraint.  These  forces  are  all  assumed  to  pass  through 
segment  mass  centers,  with  the  exception  of  the  floor  forces  and 
the  force  of  the  seat  cushion  on  segment  1.  The  floor  forces 
act  on  the  ends  of  leg  segments,  and  the  force  of  the  seat  cush- 
ion on  segment  1 passes  through  a point  midway  between  the  hip 
joints.  This  assumption  effects  a considerable  simplification 
in  the  model,  for  if  the  position  of  the  point  of  force  applica- 
tion, rp  , is  the  position  of  a segment  mass  center,  it  is  avail- 
able directly  from  the  solution  of  the  equations  of  motion. 


FORCE 

: IDENTIFICATION 

FBI, 

FB3: 

FB2 : 

BACK 

CUSHION 

HEAD  REST 

FS1 , 

FS2 , 

FS3 : 

SEAT 

CUSHION 

FF1, 

FF2: 

FLOOR 

FA1, 

FA2, 

FA3 : 

AIR  BAG 
(OPTIONAL) 

Figure  9.  External  Forces  of  Cushions,  Floor,  and  Air  Bag. 
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All  contact  forces  are  calculated  by  first  determining  the 
penetration  of  a contact  surface  on  the  occupant  into  a surface 
with  known  force-deflection  characteristics.  Using  the  seat 
cushion  force  as  an  example,  the  pertinent  dimensions  of  the  seat 
and  the  parameters  required  to  determine  the  penetration  of  the 
abdomino-pelvic  segment  {segment  1)  into  the  cushion  are  illus- 
trated in  figure  10.  Xp  and  Zp  are  coordinates  of  the  center  of 
the  contact  surface  of  segment  1,  and  is  the  radius  of  the  con- 
tact surface  in  the  (x^  - z^)  plane.  (Although  this  contact  sur- 
face is  an  ellipsoid,  cross-sections  parallel  to  the  (x^  - z^) 
plane  are  circular.  The  dimensions  of  the  contact  surfaces  will 
be  discussed  in  a later  section.)  The  position  of  the  seat  pan  is 
defined  by  its  height  Zg  above  the  origin  of  the  aircraft  co- 
ordinate system  and  the  angle  8s  that  it  makes  with  the  aircraft 
(XA  - Yft)  plane.  The  unloaded  thickness  of  the  seat  cushion  is 
t , and  the  loaded  thickness  beneath  segment  1 is  t.  Summing  the 
dimensions  in  the  Zft  direction  gives 

Zp  = Zg  + (Rx  + t)/cos  0X  + Xp  tan  6g  (36) 


Figure  10.  Seat  Cushion  Deflection. 
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Solving  equation  (36)  for  the  cushion  thickness, 

t = (Zp  ~ V COS  6g  - Rx  - Xp  sin  es 

The  deflection  of  the  seat  cushion  is  then 


(37) 


6 = t - t 


(38) 


and  the  force,  which  is  assumed  to  acr  normal  to  the  plane  of  the 
seat  pan  and  pass  through  the  center  of  curvature  of  the  contact 
surface,  is  calculated  from  an  input  curve  of  force-versus- 
def lection.  The  other  contact  forces  are  calculated  in  the  same 
manner.  A small  difference  exists  in  the  calculation  of  the  air 
bag  force  because  the  air  bag  surface  is  cylindrical,  but  the 
technique  is  basically  the  same. 

2. 1.3. 2 Restraint  System  Forces 

The  method  used  in  calculating  the  forces  exerted  on  the  body 
by  the  restraint  system  differs  considerably  from  that  described 
in  the  preceding  section  for  the  contact  forces.  The  primary  rea- 
son for  this  difference  is  that  the  restraint  forces  do  not  act  at 
any  fixed  points  on  the  occupant,  but,  rather,  the  points  of 
application  vary  with  the  restraint  system  geometry. 

Although  other  configurations  can  be  selected  by  the  user, 
the  basic  restraint  system  consists  of  a lap  belt  and  diagonal 
shoulder  strap.  The  restraint  loads  are  transmitted  to  the  occu- 
pant model  through  ellipsoidal  surfaces  fixed  to  the  upper  and 
lower  torso  segments.  These  surfaces  are  shown  in  figure  11. 

The  locations  of  the  anchor  points  A^ , A^ , and  A^  and  the  buckle 
connection  B are  determined  by  user  input  along  with  the  webbing 
properties . 

The  ellipsoidal  surfaces  are  described  by 

2 . 2 , 2 ,,  2 , 2.2  i 

xi  /ai  + n /bi  + n /ci  = 1 os) 
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and 


for  the  upper  torso. 


where 


(40) 


a0  = R„ 


c2  - V2 

and  these  body  dimensions  are  defined  in  section  2.1.4. 

The  restraint  forces  are  determined  in  the  same  runner  tor 
both  the  upper  and  1 ovci  torso.  First,  the  belt  loads  are  calcu- 
lated from  the  displacements  of  the  torso  segments,  and  the  resul- 
tant force  on  each  segment  is  then  applied  at  the  point  along  the 
arc  of  contact  between  the  belt  and  the  ellipsoidal  surface  where 
the  force  is  normal  to  the  surface. 

Explaining  this  procedure  in  further  detail  for  the  restraint 
system  configuration  shown  in  figure  11,  for  any  position  of  the 
occupant,  the  coordinates  of  the  left  shoulder,  the  hips,  and  the 
buckle  connection  R are  calculated  in  the  aircraft  reference 
frame.  The  length  of  each  of  the  three  belt  segments  (right  lap 
belt,  left  lap  belt,  and  shoulder  strap)  is  equal  to  the  sum  of 
the  free  length,  or  the  distance  between  the  appropriate  anchor 
point  and  a reference  point  on  the  hip  or  shoulder,  and  the  dis- 
tance along  the  arc  rrnm  that  hip  or  shoulder  to  the  buckle.  If 
that  length  should  exceed  the  equilibrium  (zero  load)  length 
calculated  initially  then  there  is  some  tensile  force  in  the 
belt.  The  resultant  force  on  each  segment  is  the  vector  sum  of 
the  belt  forces.  Friction  he tween  the  snoulder  belt  and  chest 
along  the  length  of  the  belt  is  taken  into  account  by  reducing 
the  load  in  the  belt  between  the  chest  and  buckle  by  a constant 
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fraction  of  the  load  in  the  free  length  between  the  belt  and  the 
body  surface.  The  resultant  force  on  the  lower  or  upper  torso 
segment  may  be  written  generally  as 


f = fx  l + Fy  j + rz  k 


(41) 


where  Fx,  F , and  Fz  are  components  in  the  local,  segment-fixed 
coordinate  system. 

To  find  the  point  on  the  segment  where  F is  normal  to  the 
surface,  consider  first  the  equation  of  an  ellipsoid: 


2.2,  2 ..2  , 2-2  , 

x /a  4-  y /b  + z /c  =1 

which  may  also  be  expressed  in  functional  form  as 


(42) 


2 2 2 2 2 2 
f (x,  y,  z)  = x/a  + y /n  + z/c  ~ 1 


(43) 


where  the  ellipsoid  can  be  regarded  as  the  level  surface  f=0  of 
the  function.  At  any  point  (x,  y,  z)  on  the  surface,  the  gradient 
of  f is  normal  tr  the  surface.  The  gradient  is  given  by 

grad  f = (2x/a2)  l + (2y/b2)  j + (2z/c2)  k (44) 

and  at  the  point  of  application  of  the  resultant  force,  grad  f is 
coliinear  with  F.  Making  use  of  the  proportionality  between  the 
components  of  the  two  vectors, 


Fv  - Cx/a2 
Fv  = Cy/b2 

u 

Fz  = Cz/c2 


where  C is  an  arbitrary  constant.  Solving  equation  (45)  for  the 
coordinates  (x,  y,  z)  and  substituting  into  equation  (42) 


ttw'iaaaia^^ 


. ,.2  0 2 0 2 
M W ?*(¥-)  > 


(46) 


which  leads  to 


(¥)’ . (¥)’  • (¥)’  • 


C2  = Fx2a2  + Fy2b2  + Fz2c2 


C = ± Vf  2a2  + F 2b2  + F 2c2 

' v \r  7 


the  point  of  application  of  F is  then 


x = FxaVC 


y = FybVC 


(47) 


(48) 


z = F c /C 
z 


with 


C 


■V3 


= -VFx2a2  + Fy2b2  + Fz2c2 


(49) 


The  negative  sign  on  C can  be  explained  by  the  fact  that  each 
coordinate  in  the  local  system  is  opposite  in  sign  to  the  cor- 
responding component  of  the  resultant  force,  or 

x > 0 if  F <0 
x 


y > 0 if  F <0 

J y 


z > 0 if  F„  < 0 


(50) 
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The  capability  of  the  belts  to  slide  relative  to  the  torso  sur- 
faces allows  simulation  of  "submarining"  under  the  lap  belt,  as 
illustrated  in  the  sequence  of  figure  12. 


Figure  12.  Occupant  Submarining. 

2.1.4  Occupant  Physical  Properties 
2. 1.4.1  Body  Segment  Dimensions 

The  basic  dimensions  of  the  occupant  segments  that  are  re- 
quired in  writing  the  equations  of  motion  are  illustrated  in 
figure  13.  The  lengths  of  the  segments  are,  in  most  cases,  effec- 
tive "link  lengths"  between  join.:  centers,  rather  than  standard 
anthropometric  dimensions  based  on  external  measurements.  These 
lengths  are  calculated  as  fixed  fractions  of  total  body  stature, 
obtained  from  references  4 and  5.  These  fractions,  along  with 
the  actual  segment  lengths  for  a 50th  percentile  male  (69.1-inc'n 
stature),  are  presented  in  table  1. 

The  distance  of  the  mass  center  of  segment  n from  the  end 
nearest  the  body  reference  point  (C^)  is  p . The  distance  be- 
tween the  mass  center  and  the  far  end  is  given  by 


31 


TABLE  1.  BODY  SEGMENT  LENGTHS 


Segment 

Symbol 

Fraction 

of 

Stature 

Length  for 
50th 

Percentile 
Male 
(in. ) 

Lower  Torso 

L1 

• 0.1714 

11.84 

Upper  Torso 

L2 

0.1571 

10.86 

Head  and  Neck 

L3 

0.1546 

10.68 

Upper  Arm 

L4 ' L6 

0.1735 

11.99 

Forearm  and  Hand 

L5 ' L7 

0.1714 

13.23 

Thigh 

L8 ' L10 

0.2399 

16.58 

Leg  and  Foot 

V L11 

0.2505 

17.31 

1/2  Hip  Breadth 

lh 

0.0491 

3.39 

Shoulder  Link 

Ls 

0.0917 

6.34 

Head  Mass  Center  Link 

I Lc 

0.0108 

0.75 

The  center  of  mass  locations,  also  obtained  from  data  of  refer 
ences  4 and  5,  are  calculated  as  fixed  fractions  of  segment 
lengths.  These  fractions  and  the  actual  dimensions  for  a 50th 
percentile  male  are  presented  in  table  2. 


TABLE  2.  SEGMENT  CENTER  OF  MASS  LOCATIONS 


Segment 

Symbol 

Fraction 
of  Segment 
Length 

Dimension 
for  50th 
"ercentile 
Male 
(in.) 

Lower  Torso 

0.4515 

5.35 

Upper  Torso 

P2 

0.4654 

5.05 

Head  and  Neck 

p3 

0.5670 

6.06 

Upper  Arm 

•V  °6 

0.4374 

5.24 

Forearm  and  Hand 

PS'  °7 

0.6770 

8.95 

Thigh 

°8'  P10 

0.4278 

7.09 

Leg  and  Foot 

°9‘  pn 

0.4264 

7.38 

2. 1.4. 2 Body  Segment  Mass 


The  masses  of  body  segments  are  calculated  in  the  same  man- 
ner • s the  lengths,  using  fractions  of  total  body  mass,  based  on 
data  of  references  4 and  5.  These  fractions  and  segment  weights 
for  a 50th  percentile  male  (161.5-pound  weight)  are  presented  in 
table  3 . 


TABLE  3.  OCCUPANT  SEGMENT  INERTIAL  PROPERTIES 

Segment 

Fraction 
of  Total 
Body  Mass 
(Mn/M) 

Weight  of 
Segment 
for  50th 
Percentile 
Male* 
(lb) 

Moments  of  Inertia 
for  50th  Percentile 
Male 

(lb-sec2-in. ) 

Xx  Iy  x* 

Lower  Torso 

0.2778 

44.86 

8.703  4.331  8.703 

Upper  Torso 

0.2264 

36.56 

3.357  2.623  2.623 

Head  and  Neck 

0.0792 

12.79 

0.311  0.311  0.201 

Upper  Arm 

0.0264 

4.26 

0.164  0.164  0.0241 

Forearm  and  Hand 

0.0214 

3.46 

0.0241  0.218  0.218 

Thigh 

0.1001 

16.17 

0.307  1.270  1.270 

Leg  and  Foot 

0.0604 



9.75 

1.192  1.192  0.120 

••Weight,  rather  than  mass  is  tabulated  here,  as  basic  units  are 
generally  more  meaningful  to  the  reader. 

2. 1.4. 3 Segment  Moments  of  Inertia 

Mass  moments  of  inertia,  computed  from  measurements  made  on 
eight  human  cadavers,  were  reported  by  Dempster  in  reference  6. 
These  moments  of  inertia  were  all  measured  with  respect  to  trans- 
verse body  axes  at  a convenient  point  of  suspension  for  each  seg- 
ment and,  subsequently,  transferred  to  the  segment  mass  centers. 
Therefore,  from  these  results  the  values  of  I can  be  extracted, 
with  the  exception  of  the  upper  torso  segment ,nwhich  will  be  dis- 
cussed further  below. 


The  cadavers  used  in  Dempster's  experiments  were  generally 
of  slight  build  possessing  an  average  weight  and  stature  consider- 
ably lower  than  the  50th  percentile.  Using  the  anthropometric 
data  on  the  United  States  civilian  population  contained  in  refer- 
ence 7 , the  weight  and  stature  for  the  50th  percentile  male  are 
161.5  pounds  and  69.1  inches,  respectively.  For  use  in  the  mathe- 
matical model,  all  of  Dempster's  values  were  adjusted  to  account 
for  total  body  size  by  multiplication  by  a factor  that  is  related 
to  the  units  of  mass  moment  of  inertia.  In  other  words,  the  ad- 
justed values  for  the  50th  percentile  male  are  calculated  accord- 
ing to 


*50  “ JD  (161‘5)  (69-1)2/KDSD2  (52) 

where  Ip  is  the  segment  moment  of  inertia  calculated  by  Dempster 
for  a population  of  average  weight  Wp  and  stature  Sp. 

Turning  to  the  upper  torso,  Dempster  subdivided  this  segment 
into  the  thorax  and  two  shoulder  segments  and  reported  the  three 
moments  of  inertia  separately.  For  use  in  the  model  a composite 
value  of  moment  of  inertia  was  calculated  as  described  below. 

Figure  14(a)  shows  the  location  of  the  center  of  mass  for 
the  thorax  anterior  to  the  vertebral  body  T-9.  Figure  14(b)  shows 
the  location  of  the  centers  of  mass  of  the  shoulder  segments  below 
rib  3.  The  lateral  location  of  the  shoulder  center  of  mass  is  un- 
important as,  for  the  two  shoulder  segments  together,  it  will  be 
in  the  midsagittal  plane.  Using  the  superior  face  of  vertebral 
body  T-l  as  a reference,  the  distances  to  the  centers  of  mass  for 
the  thorax  and  shoulders  indicated  in  figure  14  were  obtained  by 
scaling  the  dimensions  on  Dempster's  sketches.  The  length  given 
for  the  thorax  link  between  the  superior  face  of  the  centrum  of 
vertebra  T-l  and  the  inferior  face  of  the  centrum  of  vertebra  T-12 
is  15.71  percent  of  the  body  stature,  or,  for  the  50th  percentile 
male,  10.86  inches.  The  vertical  distance  to  the  center  of  mass 
of  the  thorax  is  given  as  66.1  percent  of  the  link  length  or  7.18 
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(a)  THORAX 


Figure  14.  Location  of  Centers  of  Mass  for  Upper  Torso 
Segments  (Reference  6) . 


inches  for  the  50th  percentile  male.  Using  the  spinous  process 
of  T-l,  which  is  visible  in  figure  14(b)  as  a reference,  the 
distance  from  the  superior  face  of  T-l  to  the  mass  center  of  a 
shoulder  segment  can  be  determined  to  be  approximately  28  percent 
of  the  link  length,  or  3.04  inches  for  the  50th  percentile  male. 
The  center  of  mass  for  the  combined  segment  is  then  located  5.15 
inches  or  47.4  percent  of  the  link  length  below  T-l.  The  moments 
of  inertia  can  be  combined,  using  the  parallel  axis  theorem,  as 
illustrated  in  appendix  H. 

The  moments  of  inertia  with  respect  to  the  segment  x-  and  z- 
axes  were  determined,  using  approximations  to  segment  geometry 
similar  to  the  technique  described  in  reference  8.  The  torso  and 
head  segments  were  approximated  by  ellipsoids.  Assigning  appro- 
priate anthropometric  dimensions  to  the  ellipsoid  axes,  the  ratios 

I /I  and  I /I  were  calculated  for  unit  mass.  These  ratios, 
x y z y 

n ^n  n 
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multiplied  by  the  I extracted  from  Dempster's  data,  gave  values 
of  I and  I for  the  torso  and  head  segments.  The  identical 

A £* 

n n 

procedure  was  used  for  the  extremities,  except  that  these  seg- 
ments were  approximated  by  solid  circular  cylinders.  The  de- 
tails of  the  calculations  are  presented  in  appendix  H,  and  the 
moments  of  inertia  for  a 50th  percentilfe  male  are  listed  in  table 
3. 

2. 1.4. 4 Body  Contact  Surfaces 

Twenty-three  surfaces  are  defined  on  the  body  for  calcula- 
tion of  external  forces  exerted  on  the  occupant  by  the  seat  cush- 
ions or  restraint  system  and  for  prediction  of  impact  between 
the  occupant  and  the  cockpit  interior.  These  surfaces  are  ellip- 
soids, cylinders  and  spheres,  as  shown  in  figure  15.  The  dimen- 
sions of  these  surfaces  were  obtained  from  anthropomorphic  data 
on  the  United  States  civilian  population  (reference  7). 


SURFACE  IDENTIFICATION 

1.  PELVIS 

2 . CHEST 

3 . HEAD 

4.  RIGHT  UPPER  ARM 

5.  RIGHT  FORFAPM 

6.  LEFT  UPPER  ARM 

7 . LEFT  FOREARM 

8.  RIGHT  THIGH 

9.  RIGHT  LOWER  LEG 

10.  LEFT  THIGH 

11.  LEFT  LOK1R  LEG 

12.  RIGHT  KNEE 
ii.  LEFT  KNEE 

14.  RIGHT  FOOT 

15.  LEFT  FOOT 
15.  FT  HIT  HIP 

17.  LEFT  HIP 

18.  RI JHT  SHOULDER 

19.  LEFT  SHOULDER 

20.  RIGHT  FLDCW 

21.  LEFT  ELBOW 

22.  RIGHT  HAND 

23.  LEFT  HAND 


Figure  15.  Occupant  Contact  Surfaces. 
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The  surfaces  and  the  dimensions  required  for  their  descrip- 
tion are  illustrated  in  detail  in  figure  16.  The  values  of  R^ 
for  a 50th  percentile  male  are  listed  in  table  4,  along  with 
their  fractions  of  body  stature.  The  segment  lengths  Lr  and 
center-of-mass  location  pn  were  defined  earlier  and  are  presented 
in  tables  1 and  2. 


Figure  16.  Body  Contact  Surface  Dimensions. 
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TABLE  4.  CONTACT  SURFACE  DIMENSIONS 


Surface 

Symbol 

Fraction  of 
Stature 

(Rj/S) 

Actual  Dimension 
for  50th 
Percentile  Male 
(in.) 

Pelvis 

R1 

0.0579 

4.00 

Chest 

*2 

0.0651 

4.50 

Head 

R3 

0.0579 

4.00 

Arm 

R4 ' R6 

0.0282 

1.95 

Forearm 

r5,  r7 

0.0268 

1.85 

Thigh 

V R10 

0.0514 

3.55 

Leg 

R9 ' R11 

0.0333 

2.30 

Knee 

R12 ' R13 

0.0333 

2.30 

Foot 

R14 ' R15 

0.0250 

1.73 

Hip 

R16 ' R17 

0.0515 

3.56 

Shoulder 

R18'  R19 

0.0378 

2.61 

Elbow 

R20'  R21 

0.0268 

1.85 

Hand 

R22 ' R23 

0.0297 

2.05 

I 2.. 1.4. 5 Joint  Rotation 

( 

! The  results  of  several  studies  on  the  limits  of  human  joint 

i 

met: on  have  been  published.  Two  of  these  studies  in  particular 

f 

; ' were  examined  for  applicability  to  the  occupant  model.  First  of 

l all,  Dempster's  (reference  6)  data  on  link  lengths  and  inertial 

! • properties  were  used,  as  discussed  in  preceding  sections,  so  it 

t was  considered  appropriate  to  include  his  joint  data  here.  Glan- 

ville  and  Kreezer  (reference  9)  presented  limits  of  joint  motion 
for  both  voluntary  and  forced  rotation;  their  results  appear, 
along  with  Dempster's,  in  table  5.  Definitions  of  the  various 
joint  motions  are  illustrated  in  figure  17.  Also  included  in 
table  5 are  the  rotations  recommended  for  anthropomorphic  dummy 
joints  by  SAE  Recommended  Practice  J963. 
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TABLE  5 . RANGE 

OF  JOINT  ROTATION 



Measured  Rotation  - Deg 

Human 

Body 

Component 

Glanville<9) 

and  Kreezer 

Dummy 

Dempster ^ 

Motion 

Symbol 

Motion  Description 

Voluntary 

Forced 

(SAE  J963) 

Head  - With 

A 

Dorsiflextion 

61 

77 

_ 

mm 

Respect  to 
Torso 

B 

Ventriflexion 

60  ' 

76 

- 

■98 

C 

Lateral  Flexion 

41 

63 

- 

40 

D 

Rotation 

78 

63 

- 

70 

Upper  Arm  - 

E 

Abduction  (Coronal  Plane) 

130 

137 

134 

mm 

At  Shoulder 

F 

Flexion 

180 

185 

18e 

G 

Hyperextension 

58 

69 

61 

60 

Forearm  - At 
Elbow 

H 

Flexion 

141 

146 

142 

135 

Thigh  - At 

I 

Flexion 

102 

112 

113 

120 

Hip 

J 

Hyperextension 

45 

54 

- 

45 

K 

Medial  Rotation 

- 

39 

50 

L 

Lateial  Rotation 

- 

- 

34 

50 

M 

Adduction 

- 

- 

31* 

10 

N 

Abduction 

71 

79 

53* 

50 

Lower  Leg  - 
At  Knee 

P 

Flexion 

125 

138 

125 

135 

Long  Axis  of 

SI 

Flexion 

- 

- 

- 

40 

Torso 

■9 

Hyperextension 

- 

- 

- 

30 

s 

Lateral  Flexion 

- 

- 

- 

35 

T 

Rotation 

- 

- 

- 

35 

•Transverse  plane. 

A.\l  of  the  rotations  possible  in  the  mathematical  model  are 
included  in  table  5 and  figure  17,  but  some  are,  naturally,  more 
important  than  others  in  determining  permissible  ranges  of  motion 
for  the  model.  For  the  head,  ventrif lexion  (B)  is  certainly  the 
most  important  component  of  motion  for  frontal  impact.  Dorsi- 
f lexion  (A)  may  also  be  important  for  frontal  impact,  but  the 
angles  reported  are  sufficiently  close  to  those  for  ventriflexion 
to  be  considered  the  same.  Lateral  flexion  (C)  is  certainly  less 
important  since  a pure  lateral  impact  of  an  aircraft  would  be 
rare  indeed,  and  rotation  (D)  will  have  an  insignificant  effect 
on  model  response.  Therefore,  the  limiting  rotation  0g^  (see 

section  2. 1.1. 2)  for  the  neck  joint  (i  = 2)  has  been  taken  as  the 
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limit  for  voluntary  ventri flexion,  or  60  degrees.  The  &ngle  can 
be  found  in  table  6,  along  with  the  limiting  angles  for  the 
other  body  joints.  For  all  of  the  other  angles,  flexion  is  the 
most  significant  component  for  the  type  of  motion  that  can 
usually  be  expected  to  take  place  in  a crash  environment.  There- 
fore, the  limiting  angles  were  all  taken  as  the  limits  for  volun- 
tary flexion.  Note  that,  for  the  hip  joint,  the  reference  posi- 
tion of  the  body  used  in  the  mathematical  model  includes  90-degree 
flexion.  Therefore,  this  amount  has  been  subtracted  from  the 
angle  reported  in  table  5,  which  is  defined  relative  to  the  stan- 
dard anatomical  reference  position.  Since  the  seated  position 
appears  to  aid  in  flexion  of  the  hip  joint,  the  largest  angle  in 
the  table,  the  one  given  by  SAE  J963,  was  used  in  determining 
6 which  is  thus  given  by  Bc  = 120°  - 90°  = 30°. 


TABLE 

6.  JOINT  LIMITING  ANGLES 

Angle  - & 

Si 

Joint 

Location 

(deg) 

1 

Back 

40 

2 

Neck 

60 

3,  5 

Shoulder 

180 

4,  6 

Elbow 

142 

7,  9 

Hip 

30 

8,  10 

Knee 

125 

2 . 2 SEAT  MODEL 

In  order  to  make  maximum  use  of  commonality  and  to  minimize 
input  complexity,  the  seat  model  is  divided  into  two  major  com- 
ponents. The  first,  composed  of  the  seat  pan  and  back,  retains 
the  same  configuration  regardless  of  the  particular  type  of  seat 
being  analyzed.  The  supporting  structure,  on  the  other  hand,  can 
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vary  widely  in  its  overall  geometry  and  number  of  elements.  The 
characteristics  of  each  of  the  components  of  the  seat  model  are 
discussed  below  with  attention  being  focused  first  on  the  common 
part  of  the  structure,  which  is  actually  made  up  of  three  types 
of  elements,  as  shown  in  figure  18. 


Figure  18.  Seat  Model  Components. 

The  seat  pan  is  composed  of  triangular  membrane  elements. 

The  seat  cushion  load  is  distributed  parabolically  over  the  inner 
nine  nodes.  The  sample  input  data  include  fictitious  mechanical 
properties  for  these  elements  that  have  been  determined  to  best 
simulate  a fabric  or  spring-supported  seat  pan.  A sheet-metal 
seat  pan  would  be  represented  by  using  actual  material  properties 
as  input  data. 


43 


A seat  pan  frame  formed  of  16  beam  elements  transmits  the 
seat  pan  loads  to  the  supporting  structure.  Also  applied  to  this 
frame  are  loads  from  the  seat  back,  the  occupant's  thighs,  and 
the  lap  belt,  should  the  user  wish  the  belt  attached  to  the  seat. 

The  seat  back  is  made  up  of  three  simple  beams  of  '■•niform 
cross  section.  The  back  cushion  loads  are  assumed  to  be  distri- 
buted along  the  sides  cf.  cne  seat  back  and  are  thus  transmitted 
to  the  seat  pan  frame.  Should  the  shoulder  harness  be  attached 
to  the  seat,  its  load  is  applied  at  the  midpoint  of  the  transverse 
beam  that  forms  the  top  of  the  seat  back.  Bending  of  these  beams 
can  alter  the  seat  geometry  during  the  crash  event,  and  bending 
failure  at  the  connection  to  the  seat  pan  frame  can  be  predicted, 
based  on  an  input  value  of  ultimate  stress  for  the  frame  material. 

2.2.1  Seat  Analysis 

The  behavior  of  a typical  occupied  seat  in  an  impact  situa- 
tion is  characterized  by  three  distinct,  though  overlapping 
phases.  Initially  the  seat  pan  is  only  slightly  deformed  and 
comparatively  flexible.  The  resistance  of  such  a "flat"  mem- 
brane to  normal  displacement  arises  primarily  from  its  current 
state  of  stress  and  current  geometry.  Once  the  seat  pan  is  suf- 
ficiently deformed  to  be  of  comparable  stiffness  to  the  frame, 
the  second  phase  of  behavior  is  pursued.  Herein  the  framework 
and  the  pan  continue  to  deform  elastically  with  neither  contribu- 
tion necessarily  dominant.  The  flexure  of  the  supporting  frame 
acts  to  soften  the  pan  stiffness  while  continued  stretching  of 
the  pan  counterbalances  this  influence.  In  the  final  phase, 
frame  components  may  become  plasticized,  thereby  becoming  in- 
capable of  carrying  greater  load.  Plastic  hinges  are  introduced, 
causing  redistribution  of  the  loads.  Redistribution  continues 
until  sufficient  hinges  are  formed  to  allow  the  structure  to  de- 
form as  a kinematic  mechanism,  at  which  point  collapse  occurs. 

The  present  analysis  accounts  for  the  above  described  be- 
havior although  no  attempt  is  made  to  discretize  the  separate 
phases.  All  phenomena  are  considered  throughout  the  analysis, 
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thus  providing  gradual  transition  from  initiation  of  deformation 
to  termination,  possibly  by  total  collapse. 

The  data  handling  for  the  seat  analysis  is  designed  to  ef- 
ficiently utilize  limited  computer  storage  capacity  by  not  allow- 
ing large  matrices,  generated  for  various  components,  to  reside 
in  core  simultaneously.  The  general  flow,  described  in  more  de- 
tail subsequently,  involves  first  forming  the  structural  matrix 
describing  the  particular  supporting  structure.  This  matrix  is 
modified  during  formulation  to  account  for  any  plastic  hinges  or 
floor  connection  failure.  Once  formed,  the  matrix  is  copied  to 
auxiliary  storage  and  numerically  collapsed  to  the  seat  pan  frame 
intersection  nodes.  The  seat  pan  stiffness  matrix  is  then  cre- 
ated, added  to  the  leg  frame  residue,  and  the  resulting  combined 
matrix  saved  on  auxiliary  storage  until  further  collapsed  to  the 
seat  pan  intersection. 

Similarly,  the  seat  pan  matrix  is  formed  and  solved  for  the 
increments  of  displacement.  Frame,  leg,  and  bar  matrices  are  re- 
trieved and  used  to  calculate  bar  forces  and  floor  reactions. 
Updated  displacement  data  are  stored  to  define  current  position 
in  the  subsequent  call  to  seat  routine. 

Discussions  of  analysis  techniques  performed  in  the  program, 
which  are  applicable  to  one  or  more  of  the  seat  types,  are  pre- 
sented below. 

2. 2. 1.1  Bar  Stiffness 

A straight  bar  of  constant  cross  section,  as  shown  in  figure 
19,  is  taken  as  the  first  element.  For  small  elastic  deformation 
the  customary  approximations  used  in  finite-element  developments 
are  unnecessary  becaus°  the  force-displacement  relationships  are 
well  known.  The  six  components  of  displacement  at  each  end  are 
related  to  the  corresponding  generalized  forces  by  the  matrix 
expression 


{P}  = {K}  {p} 


(53) 
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Figure  19.  Space  Frame  Element  in  Elemental 
Coordinate  System. 


where  { P > = the  12  x 1 column  vector  of  nodal  forces 

{ p } = the  12  x 1 column  vector  of  nodal  displacements 

{K}  = the  12  x 12  stiffness  matrix 

In  expanded  notation  the  equivalent  expression  is  given  by  equa- 
tion (54) 


where 


M. 


V 


A = area  of  cross  section 
E = modulus  of  elasticity 
I = moment  of  inertia 
G = modulus  of  rigidity 
J = torsion  constant 
L = length  of  member 
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0 » (0  , ) # (0  ) = rotation  about  X,  (Y) , (Z)  axis  at 

xi  yi  zi  node  i 
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2. 2. 1.2  Rigid  Body  Transformation 

In  general,  actual  seat  supporting  structure  attachments 
might  only  rarely  coincide  with  natural  integral  division  points 
of  the  seat  pan  frame.  A rigid  body  transformation  technique  is 
employed  to  accommodate  this  "mismatch",  without  increasing  the 
problem  sivse,  by  the  obvious  technique  of  adding  nodes.  By  this 
method  nodal  displacements  and  bar  end  displacements  are  related 
through  kinematics,  resulting  in  a local  approximation  of 
stresses,  but  providing  full  equilibrium  and  continuity. 


The 


Figure  20  illustrates  the  geometry  and  nomenclature, 
bar  end  B is  eccentric  from  the  node  N by  the  distances 

(ex'  ey ' ez) 


Figure  20.  Geometry  of  Rigid-body  Transformation. 


The  bar  end  displacements  are  related,  in  a purely  kine- 
matic sense,  to  the  noda1  displacements  by  the  matrix  equation 


(eB)  = [A]  {eN> 


(55) 


and  {eB(N)}  represents  the  6x1  vector  of  displacement  at  B (N) 
and  [A]  is  the  rigid  body  relation  defined  by 


(56) 
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The  transpose  of  iA]  represents  the  same  rigid  body  trans- 
formation applied  to  the  corresponding  force  components. 

PB  * (57) 


Th'  transformed  stiffness  matrix  is  then  given  by  the  con 
gruent  transformation 


K'  = TT  K T (58) 

where  [T]  is  formed  from  [A] 


T12x12  " 


2. 2. 1.3  Matrix  Reduction 

In  order  to  use  available  computer  storage  economically  it 
is  necessary  to  reduce  larger  order  matrixes  for  components  to 
their  influence  on  freedoms  common  to  other  components.  The  gen- 
eral technique  explained  subsequently  performs  this  operation. 

The  subscript  i refers  to  the  variables  to  be  eliminated  while  j 
refers  to  those  retained.  The  combined  (i  and  j)  matrix  is  par- 
titioned into  the  indicated  submatrixes. 


K = 


K . . I K . . 

3.1  | 13 

I 

I 

Kji  ! 


The  corresponding  equilibrium  equations  are 


Kii  ri  + Kij  rj  = Ri 


(60) 


(61) 
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(62) 


Kij  ri  + Ki-J  • 'S 


Solution  of  equation  (61)  gives 


ri  • Ku'1  IRi  - Kij  rj) 


(63) 


Substitution  of  equation  (63)  into  equation  (62)  gives 


[K..  - K..  K.."1  K..]  r.  = R.  - K.  . K . . _1  R. 

33  3i  11  13  3 3 13  11  1 (64) 


The  bracketed  quantity  on  the  left  side  of  the  equation  is 
the  modified  or  reduced  stiffness  matrix;  a modified  load  vector 
is  represented  on  the  right  side. 

2. 2. 1.4  Plastic  Modification 

A plastic  hinge  is  introduced  at  any  point  in  the  analysis 
when  the  bending  moment  at  a bar  end  reaches  a predetermined 
limiting  value.  This  requires  a modification  of  the  bar  stiffness 
matrix  to  remove  the  capability  to  resist  further  moment.  The 
modification  is  performed  as  follows. 

The  stiffness  matrix  presented  previously  represents  the 
self-equilibrating  set  of  forces  at  the  bar  ends  induced  by  the 
end  displacements.  The  physical  condition  that  one  of  these  com- 
ponents be  zero  is  accomplished  by  the  matrix  reduction  of  the 
previous  section,  where  in  this  application  K refers  to  the 
12x12  bar  stiffness  matrix,  i is  the  component  of  zero  moment, 
and  j the  remaining  nonplasticized  degrees  of  freedom. 

The  ii  term  is  replaced  by  a small  number  tor  convenience 
in  computation  to  prevent  an  apparent  mechanism  from  developing. 

2. 2. 1.5  Triangular  (Membrane)  Element  Stiffness 

Elastic  Stiffness  - The  triangular  panel  of  figure  21  is 
considered.  The  bending  stiffness  is  ignored  and  the  state  of 
membrane  stress  is  assumed  constant  throughout  the  element. 
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Figure  21.  Triangular  Panel  Element. 

Nodal  points  are  the  three  corner**.  The  relationship  between 
representative  nodal  force  and  in-plane  nodal  displacements  (or 
stiffness  matrix)  is 
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Geometric  Stiffness  - The  geometric  stiffness  of  a triangu- 
lar panel  is  presented  subsequently.  Physically,  it  represents 
the  forces  required  to  displace  (rotate)  a stretched  membrane 
normal  to  its  plane. 

The  stiffness  is  expressed  in  matrix  form  as 


t 


dA 


(66) 


where 


0 (a  ) 
x y 


t 

A 


w 


the  direct  stress  in  the  x(y)  direction 

the  membrane  shearing  stress 

plate  thickness 
area 

displacement  normal  to  the  plane  of  the  plate 


and  x,  y are  the  in-plane  position  ordinates;  w , w are  slopes 

x y 

in  the  subscripted  directions.  With  the  assumption  that  the 
average  slopes  are  sufficiently  representative  of  the  detailed 
deformation,  the  derivatives  w^,  w^  are  expressed  directly  in 
terms  of  the  corner  displacements  normal  to  the  plate  surface. 


2. 2. 1.6  Coordinate  Transformations 


The  elemental  stiffness  matrix  formulations  are  most  con- 
veniently carried  out  in  the  natural  local  coordinates  of  the 
element. 


The  congruent  rigid  body  transformation,  in  reality,  is  only 
a rather  special  case  of  generalized  coordinate  transformation. 

In  the  general  case,  wherein  T represents  the  relation  between 
coordinate  systems,  the  congruent  transform  becomes 


K 


T 

T K T 


(67) 
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For  the  particular  case  of  axis  rotation 


and 


T - 


(68) 


t = 


m 


in 


m_ 


n 


n. 


(69) 


where  l , i , Z are  the  direction  cosine  of  global  X in  the 
y the  local  xyz  system 

mx,  my,  mz  are  the  direction  cosines  of  global  Y 

n , n , are  the  direction  cosines  of  global  Z 
x y - 

2. 2. 1.7  Master  Matrix  Formulation 

Elemental  matrices  are  combined  by  the  direct  stiffness 
method  into  a "master"  matrix  for  the  seat  structure.  This  pro- 
cedure can  be  represented  by  the  matrix  equation 


K = [B]1  [Keltl  [B] 


(70) 


where  Ke^t  contains  the  element  stiffness  matrices  as  submatrices 
on  the  main  diagonal  and  zeros  elsewhere. 

k. 


K 


elt 


(71) 
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and  B is  a Boolean  transformation  matrix  relating  the  local  and 
master  degrees  of  freedom. 

Computationally,  this  operation  never  actually  materializes 
because  it  would  physically  amount  to  adding  the  local  stiffness 
matrices  into  the  global  system  in  proper  order.  This  is  best 
accomplished  one  element  at  a time  by  direct  addition. 

2.2.2  Seat  Types 

The  configuration  of  the  supporting  structure  can  be  varied 
in  order  to  describe  a particular  seat.  The  only  information 
required  for  input  involves  seat  dimensions  and  material  proper- 
ties. The  simplicity  of  the  input  data  is  possible  because  the 
geometry  of  four  different  types  of  seat  configurations  is  in- 
cluded in  the  program.  These  four  seat  configurations,  which  are 
capable  of  simulating  nearly  every  seat  found  in  light  aircraft, 
are  defined  below  and  illustrated  in  figure  22. 


Figure  22.  Seat  Types. 
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Type  1.  The  first  seat  type  is  the  most  general  and  is  re- 
presentative of  the  greatest  number  of  seat  configurations. 
An  example  is  shown  in  figure  22 (a) . The  supporting  struc- 
ture which  need  not  be  symmetric,  consists  of  an  arbitrary 
number  of  beam  elements  that  connect  up  to  six  nodes  on  the 
aircraft  floor  with  any  of  the  10  nodes  on  the  outer  edges 
of  the  seat  pan  frame,  as  well  as  with  each  otnor.  As  an 
example,  the  type  1 seat  shown  in  figure  18  has  a support- 
ing structure  consi^.  ;ing  of  six  bars  that  connect  four  floor 
nodes  with  four  frame  nodes  and  with  each  other. 

The  primary  restriction  on  the  response  of  the  type  1 seat 
is  that  all  frame  nodes  are  fixed  to  resist  rotation,  as  in 
the  case  of  welded  joints.  However,  at  any  point  in  the 
analysis,  should  the  bending  moment  at  a bar  end  reach  a 
predetermined  limiting  value,  a plastic  hinge  is  introduced. 
The  limiting  moments  are  based  on  input  values  of  yield 
stress  and  provide  a means  for  simulating  bar  failure. 

Type  2.  The  second  seat  type  has  up  to  14  beam  elements  in 
the  supporting  structure.  Bars  may  connect  several  pinned- 
type  joints  on  the  seat  legs,  but  the  number  of  floor  nodes 
is  restricted  to  four.  An  example  is  shown  in  figure  22(b). 

Type  3.  This  seat-  configuration  has  a substantially  more 
complex  structure  than  the  previous  two.  As  shown  in  figure 
22(c),  the  seat  is  vertically  adjustable  with  adjustment 
provided  by  motion  of  the  seat-supporting  frame  along  rela- 
tively rigid  guide  tubes. 

Type  4.  This  seat  is  another  vertically  adjustable  config- 
uration with  adjustment  provided  by  rotation  of  crank  mecha- 
nisms, as  shown  in  figure  22(d). 
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3.0  SIMULATION  COMPUTER  PROGRAM 

The  digital  computer  program  based  on  the  occupant  and  seat 
models  described  in  section  2.0  is  called  Seat  Occupant  Model- 
Light  Aircraft  (SOM-LA) . It  has  been  written  entirely  in 
FORTRAN  IV  to  ensure  a high  degree  of  compatibility  with  various 
digital  computer  systems.  During  development  the  program  has 
been  run  on  both  UNIVAC  1108  and  CDC  6600  computer  systems. 

The  elements  of  the  program  can  be  considered  in  terms  of 
three  general  operations: 

• Input  and  Initialization 

• Solution 

• Output 

which  are  summarized  below  and  discussed  in  detail  in  the  sections 
following.  The  general  flow  of  the  program  is  illustrated  in 
figure  23.  Input  data  describing  the  occupant  and  crash  condi- 
tions are  read  first.  If  the  user  requests  output  of  the  pre- 
diction of  impact  between  the  occupant  and  the  aircraft  interior, 
the  coordinates  defining  the  cockpit  surfaces  are  read.  Finally, 
the  seat  data,  either  simple  dimensions  describing  a rigid  seat 
model  or  detailed  design  data  on  a flexible  seat,  are  provided. 
Based  on  the  input  data,  the  values  of  constants,  such  as  occu- 
pant dimensions  and  properties  are  calculated,  and  the  initial 
position  of  the  occupant  is  determined. 

The  solution  loop  is  entered  for  the  first  time  with  the  air- 
craft initial  velocity  and  the  occupant  initial  position.  At 
each  subsequent  entrance  to  the  loop  the  current  aircraft  dis- 
placement, velocity,  and  acceleration  components  are  calculated. 
The  equations  of  motion  for  the  occupant  are  set  up  and  solved. 

If  a flexible  seat  is  being  used,  the  forces  applied  to  the  seat, 
such  as  the  cushion  forces,  are  provided  to  the  seat  routines  for 
computation  of  seat  displacements.  At  time  increments  equal  to  a 
predetermined  print  interval,  the  output  variables  requested  by 
the  user  are  stored  for  printing  after  completion  of  the  solution. 
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3.1  PROGRAM  INPUT  AND  INITIALIZATION 

Input  data  are  read  by  the  program  in  the  following  seven 


blocks : 

1. 

Simulation  control  information 

2. 

Occupant  description 

3. 

Restraint  system  description 

4 . 

Cushion  properties 

5. 

Cockpit  description 

6. 

Crash  conditions 

7. 

Seat  design  information 

3.1.1  Simulation  Control  Information 

The  first  block  of  data  contains  the  information  required 
for  controlling  execution  of  the  program.  The  initial  time  step 
for  integration  of  the  equations  of  motion,  the  total  length  of 
the  simulation,  the  number  of  cases  to  be  run,  the  system  of 
units  (SI  or  English) , and  identification  of  the  desired  output 
are  provided  here. 
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3.1.2  Occupant  Description 

Because  it  has  been  assumed  that  the  principal  user  of  this 
program  is  interested  primarily  in  the  seat  or  restraint  system, 
a minimum  of  information  is  required  to  describe  the  occupant. 

Data  include  a parameter  that  defines  the  occupant  type  - human 
or  dummy  - and  the  size  of  that  occupant.  For  a human  occupant 
the  stature  and  total  weight  are  required;  for  a dummy,  the  per- 
centile (95th  male,  50th  male,  or  5th  female) . 

3. 1.2.1  Occupant  Properties 

The  dimensions  and  inertial  properties  of  the  11  body  seg- 
ments are  determined  within  the  program,  as  discussed  in  section 
2.1.4. 

3. 1.2. 2 Occupant  Initial  Position 

The  initial  position  of  the  aircraft  occupant  is  computed 

from  the  input  parameters  shown  in  figure  24.  It  is  assumed 

that  the  occupant  is  seated  symmetrically  with  respect  to  the 

aircraft  - ZA  plane  or,  equivalently,  that  the  seqment  fixed 

y -axes  are  all  parallel  to  the  Y-axis.  The  angular  coordinates 
n.  a 

(i  = 1,  2,  3,  4)  define  the  rotation  of  segments  1-4  relative 
to  the  Z axis  and,  because  of  the  symmetry  condition,  segment  6 

A 

is  parallel  to  segment  4.  The  angle  a describes  the  position  of 
the  forearms  relative  to  the  upper  arms,  and  is  the  initial  value 
of  ac  and  a_.  The  distance  X„  is  the  initial  X-coordinate  of  the 

Of  fl 

heels  (the  inferior  ends  of  segments  9 and  11) . The  procedure 
described  below  consists  of  seating  the  occupant  in  such  a posi- 
tion that  static  equilibrium  is  achieved  among  the  forces  exerted 
by  the  seat  cushion,  floor,  and  either  the  restraint  system  or 
the  back  cushion. 

The  first  step  in  determining  the  initial  position  involves 
calculating  the  Euler  angles  for  the  torso,  head,  and  arm  segments, 
since  this  procedure  does  not  require  consideration  of  the  forces 
due  to  the  cushions  and  floor.  Because  the  input  parameters 
illustrated  in  figure  24  define  the  position  of  the  occupant  in 
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a 


Figure  24.  Initial  Position  Input  Parameters. 


| 1 the  aircraft  coordinate  system,  the  orientation  of  the  aircraft 

must  be  described  in  the  inertial  system.  For  an  aircraft  in 
level  flight  with  zero  pitch,  roll,  and  yaw  it  is  assumed  that 
| the  aircraft  coordinate  axes  <XA,  Yft,  Zft)  are  parallel  to  the 

l fixed  coordinate  axes  (X,  Y,  Z)  at  the  initial  time  (t  = 0) . A 
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general  orientation  of  the  aircraft  reference  frame  is  obtained 
by  the  same  sequence  of  rotations  defined  in  section  2. 1.1.1  for 
the  body  segments.  Defining  the  rotations 

to, : Yaw 

A 

0.:  Pitch 

A 

: Roll 

A 

the  orientation  of  the  aircraft  relative  to  the  inertial  system 
is  described  by  the  coordinate  transformation 


X X. 

A 

Y * A Ya 

Z ZA 

where  the  elements  of  [A]  are 
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(73) 
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The  rotation  of  body  segment  n relative  to  the  aircraft,  re- 
membering that  the  symmetry  condition  requires  that  yn  is  parallel 
to  Ya,  is  described  by 


cos  y 0 sin  y 
n n 


-sin  y 0 cos  y 
n n 


(74) 


Combining  equations  (72)  and  (74)  results  in  the  angular  relation- 
ship between  the  local  coordinate  system  of  segment  n and  the 
inertial  system  expressed  by  the  following  transformation,  which 
is  a function  of  the  input  Yn  and  the  aircraft  pitch,  roll,  and 
yaw: 


v 


where  [Bn]  is  given  by 


M - 

A 

cos  y 
'n 

0 

0 sin  y„ 

n 

1 0 

_ « 

. -sin  yn 

0 cos  Yn 

(75) 


(76) 


so  that  its  elements  are 


B 
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= A11 
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B13n  = A11  sin  Yn  + A13  008  Yn 

B2l"  * A21  005  Yn  - s23  sin  Yn 

n ® — a 
“22  A22 

B23n  ' A21  sin  Yn  + A23  C0S  Yn 

B3l"  " A31  oos  Yn  ’ A33  sin  Yn 

B n - R 

“32  A32 

B33n  ’ A31  sin  Yn  + A33  cos  Yn  (77) 

Comparison  of  equation  (75)  with  equation  (1)  points  out  that  the 
transformation  matrices  [Tn]  and  [Bn]  are  equivalent.  Because 
[Tn]  is  a function  of  the  Euler  angles  for  segment  n,  equating 
the  elements  of  [Tn]  and  [Bn]  through 

[Tn]  = [Bn]  n = 1,  2,  3,  4,  6 , 8 , 10  (78) 

permits  calculation  of  the  initial  values  of  the  generalized  co- 
ordinates from  input  parameters  ya  and  \pA,  0 A,  and  <J>A*  The  pro- 
cedure as  used  in  Program  SOM-LA  is  outlined  below. 

First  0 is  determined  as  follows: 
n 

T31  = B31  or  'sin  0n  = B31 

gives 


9n  - "i"'1 


The  cosine  is  then  found  by 


— n 

cos  0n  = cos  [sin"'L  (-B3P ) 


(79) 
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so  that  ipn  can  be  determined: 

T11  = B11  or  008  008 

gives 

••1  Y\ 

= cos”  (B^/cos 
and,  for  determination  of  <J>n, 

T33  = B33  or  cos  6n  008 

gives 

<j>n  = cos”1  (B^/cos 

Equations  (79)  through  (81)  are  used  for  segments  1,  2,  3, 
and  4;  the  symmetry  requirement  provides  the  Euler  angles  for 
segment  6,  which  are  equal  to  those  for  segment  4.  At  this 
point  the  generalized  coordinates  q4  through  q2Q  have  been  de- 
termined. The  next  step  involves  seating  the  occupant  and  calcu- 
lating X^,  Y^,  and  Z1  (q1#  q2,  and  q^)  from  static  equilibrium. 

Because  the  problem  of  seating  the  occupant  is  statically 
indeterminate , certain  simplifying  assumptions  are  made.  The 
first  assumption,  which  is  approximately  correct  for  typical 
seating  positions,  is  that  15  percent  of  the  occupant's  weight  is 
supported  by  the  floor.  In  other  words,  85  percent  is  supported 
by  the  seat  cushion  and,  depending  on  the  aircraft  attitude,  the 
restraint  system  or  the  back  cushion. 

A first  approximation  to  the  initial  position  is  made  for  the 
assumption  of  level  flight  (@A  «=  iJ>A  = <J>A  = 0) . The  cushion  forces 
act  on  the  body  as  shown  in  figure  25,  where  it  is  assumed  that 
fifteen  percent  of  the  occupant  weight  is  supported  by  the  floor, 
as  discussed  in  the  preceding  paragraph.  Summing  forces  gives 


9n  “ BU 


V 


(80) 


= b,5 

Yn  33 


V 


(81) 
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Figure  25.  Forces  Acting  on  Occupant  Torso 
(Level  Flight) . 


Fx  : Fg  cos  9g  - Fg  sin  0g  = 0 

aa 

Fz  : Fg  sin  0g  + Fg  cos  0g  = W' 

A 

which  can  be  solved  for  the  cushion  forces: 

Fg  ~ W’  COS  0g/COS  (0g  - 0g) 

Fg  = W sir.  0g/cos  (0g  - 0g) 

Dimensional  considerations  permit  the  coordinates  of  point  P to  be 
written  as  functions  of  the  thicknesses  tg  and  tfi  of  the  compressed 
seat  and  back  cushions,  respectively. 
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(82) 


(83) 


zp  * zs  + (Rx  + tg)/cos  0g  + xp  tan  0g 

xp  = XS  + (R1  + tB>/cos  9b  " zp  tan  0b 
which  can  be  solved  for  X and  Z to  qive 

p p 

Xp  “ *fl  cos  6S  " f 2 sin  0B)/COS  (9b  " 9s* 

zp  = (fx  sin  0g  + f2  cos  0B)/cos  (0fi  - 0g) 

where  fi  * Ri  + tB  + Xg  cos  0fi 


£2  = Ri  + ts  + zs  cos  es  (8S) 

Since  the  force-deformation  characteristics  of  the  cushions  are 
known  from  input  data,  the  compressed  thicknesses  tg  and  tg  can 
be  calculated  from  equation  (83).  These  values,  when  used  in 
equation  (85) , give  the  coordinates  of  point  P for  the  first 
approximation  of  level  flight.  The  equilibrium  (zero-load)  lengths 
of  the  lap  belt  and  shoulder  belt(s)  are  calculated  for  the  body 
in  this  position. 

Next,  the  aircraft  is  rotated  to  the  attitude  specified  by 
the  input  conditions  of  pitch,  roll,  and  yaw.  Nose-up  pitch  will 
tend  to  load  the  back  cushion,  and  the  analysis  will  be  the  same 
as  that  described  abo^e  for  level  flight,  except  that  the  W' 
vector  in  figure  25  will  have  a component  in  the  XA  direction. 

Nose-down  pitch,  on  the  other  hand,  will  tend  to  load  the 
restraint  system.  An  iterative  procedure  is  used  to  determine 
the  correct  position  for  this  case.  Referring  to  figure  26,  sum- 
ming forces  gives  a set  of  transcendental  equations 


Fx  : W'  sin  0A  - Fg  sin  0g  - F^  cos  0L  = 0 

Fg  : -W'  cos  0A  + Fg  cos  0g  - Fl  sin  0L  = 0 


(86) 
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Figure  26.  Forces  Acting  on  Occupant  Torso 
(Nose-Down  Attitude)  . 

where  0L  is  the  angle  between  the  floor  and  the  projection  of  the 
lap  belt  on  the  - ZA  plane.  An  initial  estimate  to  0L  is  made 
using  the  body  position  calculated  by  equation  (85)  for  the  level 
flight  assumption.  The  angle  is  defined  according  to 


The  forces  in  the  seat  cushion  Fg  and  the  lap  belt  are  deter- 
mined using  this  value  of  0^  in  equation  (86) . From  the  input 
force-deformation  characteristics  for  the  seat  cushion  and  lap 
belt,  the  deformations  6g  and  6^  are  calculated.  Thest.  deforma- 
tions are  used  to  determine  new  values  of  Xp  and  Zp;  this  proce- 
dure amounts  to  permitting  the  body  to  further  compress  the  seat 
cushion  and  slide  forward  into  the  lap  belt.  Following  through 
the  procedure,  the  new  length  for  one  side  of  the  lap  belt  i - 

LI,  = LLe  + 6L  (88) 

where  LLe  is  the  equilibrium  length.  The  new  value  of  Xp  is  given 
by 


XP  - XL  + l(LL  ' V2  - «P  - V 2lV2  008  6L  (89) 

where  LH  is  one-half  the  hip  breadth  and  Yp  is  the  Y-coordinate  of 
the  right  hip  in  the  aircraft  system.  The  new  value  of  Zp  is  com- 
puted for  the  new  cushion  thickness  tg  using  equation  (84) , which 
is  repeated  here  for  continuity: 

zp  = Zg  + (^  + tg)/cos  0g  + Xp  tan  0g 

The  new  occupant  position,  determined  by  equations  (89)  and 

(84)  is  used  in  equation  (87)  to  recalculate  the  lap  belt  angle 

0T , and  the  procedure  is  repeated  until  two  consecutive  values  of 
h 

Xp  agree  by  less  than  5 percent.  The  coordinates  of  the  mass 
center  of  segment  1 (X^  Y^,  Z^)  are  then  calculated  from  Xpf  Yp, 
and  Zp. 

At  this  point  the  generalized  coordinates  through  q2Q 
have  been  determined.  The  final  task  will  be  to  determine  the 
coordinates  for  the  legs.  Referring  to  figure  27,  the  angles 
Y8  and  0^  can  be  found  from  simple  geometric  relationships  among 
the  dimensions  shown.  The  Euler  angles  \p g,  0g,  and  <j>g  are  ob- 
tained from  y8  by  using  equations  (79)  through  (81),  and  the 
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Figure  27.  Leg  Position. 

corresponding  coordinates  for  segment  10,  by  symmetry.  The  knee 
angles  are  given  by 


“m  = " - 6k 


m = 9 , 11 


to  complete  the  initialization  of  the  28  generalized  coordinates. 

An  alternative  procedure,  which  is  included  in  program  SOM-LA 
to  permit  re-start  from  the  final  conditions  of  a previous  run, 
bypasses  the  entire  initialization  procedure  described  here  and 
uses  input  values  of  q1  - q2Q. 

3.1.3  Restraint  System  Description 


The  restraint  system  used  in  the  simulation  may  consist  of 
a lap  belt  alone  or  combined  with  a single-  or  double-strap  shoul- 
der harness.  The  webbing  force-elongation  curve  is  approximated 
by  three  linear  segments,  which  are  described  by  input  of  four 
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points  on  the  curve.  The  force  is  computed  by  linear  interpola- 
tion in  this  table , as  described  in  section  3.2.1.  The  slack  in 
the  webbing  is  also  provided  by  input  in  units  of  length. 

The  anchor  points  for  the  lap  belt  and  shoulder  harness  are 
located  by  input  of  rectangular  coordinates  in  the  aircraft  ref- 
erence system.  For  a double-strap  shoulder  harness  the  buckle, 
or  point  of  connection  to  the  lap  belt,  is  assumed  located  on  the 
m:d-point  of  the  lap  belt.  Tor  a single  shoulder  belt,  which 
may  pass  over  either  the  left  or  the  right  shoulder,  an  input 
parameter  locates  the  buckle  by  the  length  of  webbing  between 
the  buckle  and  the  lac  belt  anchor  point.  This  length  may  be 
zero  if  the  buckle  attaches  directly  to  a rigid  anchor  point. 

3.1.4  Cushion  Description 

Input  of  a table  of  four  forces  and  deflections  describes 
the  cushion  characteristics.  The  equilibrium  (zero  load)  thick- 
ness for  both  the  seat  and  back  cushions  are  also  given. 

3.1.5  Cockpit  Description 

For  prediction  of  impact  between  the  occupant  and  the  cock- 
pit interior,  ten  plane  surfaces  are  used  to  describe  the  cockpit. 
As  shown  in  figure  28,  six  of  these  surfaces  are  normal  to  the 
- ZA  plane  and  four  are  normal  to  the  YA  - ZA  plane.  The 
first  five  planes  can  be  used  to  describe  the  environment  of  a 
crew  seat,  in  which  case  they  represent  the  firewall,  instrument 
panel,  and  windscreen,  or,  for  analysis  of  a passenger  seat,  they 
can  be  rearranged  to  describe  a seat  back.  Input  data  include  X 
and  Z coordinates  to  define  planes  1-6  and  Y and  Z coordinates 
for  planes  7-10. 

3.1.6  Crash  Conditions 

The  aircraft  crash  conditions  are  defined  by  the  initial 
velocity  and  attitude  and  the  acceleration  as  a function  of  time. 
Six  components  of  velocity  are  required:  three  translational  in 

the  aircraft  coordinate  system  (Vv  , Vv  , V_  ) and  the  yaw,  pitch, 
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Figure  28.  Planar  Surface  Approximation 
to  Aircraft  Interior. 


FRONT  VIEW 


and  roll  rates  (tjiA,  0A,  <$>A) . Each  of  the  six  acceleration  com- 
ponents, which  define  the  acceleration  of  the  aircraft  coordinate 
system,  is  described  by  sixteen  points  in  time  and  acceleration. 
An  example  of  an  approximation  to  an  actual  acceleration  pulse 
is  illustrated  in  figure  29.  Although  many  of  the  higher  fre- 
quency oscillations  observed  in  the  actual  pulse  probably  contri- 
bute little  to  the  overall  response  of  the  occupant,  the  use  of 
a large  number  of  points  reduces  the  effect  of  the  investigator's 
subjectivity  in  the  approximation. 

3.1.7  Seat  Design  Information 

The  input  data  required  to  describe  the  seat  consist  of 
dimensions,  material  properties,  and  floor  attachment  strengths, 
as  described  earlier.  For  possible  use  in  restraint  system  or 
cabin  configuration  analyses  where  seat  response  may  be  unimpor- 
tant or  seat  design  unknown,  a rigid  seat  model  can  be  selected. 
Input  data  for  the  rigid  seat  option  consist  only  of  locations 
of  the  seat  pan  and  seat  back. 
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Figure  29.  Approximation  to  Acceleration  Pulse. 

3.2  PROGRAM  SOLUTION  PROCEDURE 

The  first  operation  in  each  solution  step  includes  the  calcu- 
lation of  new  values  for  the  aircraft  acceleration  components  and 
their  subsequent  integration  to  obtain  aircraft  velocity  and 
displacement  components.  Then  the  matrix  form  of  the  equations 
of  motion,  using  equation  (19) 

[A  (q)  ] {q}  = {B (q ,q)  } + (P(q)}  + (R(q,q)}  + (Q(q,q)) 

are  set  up  for  solution  and  solved,  as  discussed  below. 
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3.2.1  Setup  of  Equations  of  Motion 

The  elements  of  [A],  {B},  {P},  and  {R}  are  calculated  from 
the  expressions  presented  in  appendixes  D through  F,  using  the 
current  values  of  the  generalized  coordinates  and  velocities. 

The  elements  of  {Q},  which  is  the  vector  of  generalized  external 
forces,  are  calculated,  as  discussed  in  section  2.1.3.  The  ex- 
ternal forces  depend  on  displacements  of  the  aircraft,  which 
determine  the  motion  of  the  seat,  floor,  and  restraint  system 
anchor  points  relative  to  the  body.  From  these  displacements 
new  deflections  of  the  cushions,  floor,  and  restraint  system  are 
calculated.  The  forces  are  then  calculated  by  linear  interpola- 
tion in  a table  of  forces  and  deformations  are  provided  as  input 
data. 

The  model  used  for  calculating  all  forces  is  illustrated  in 
figure  30,  which  shows  the  force-elongation  characteristics  for 
a typical  restraint  system  webbing.  The  experimental  curve  is 
approximated  by  three  linear  segments.  Deformation  between 
points  1 and  2 is  considered  elastic,  so  that  ux loading  would 
proceed  down  the  loading  curve.  However,  beyond  point  2 a de- 
crease in  deformation  will  cause  the  member  to  unload  along  a 
fourth-order  curve  between  the  point  where  the  unloading  starts 
and  the  origin.  If  the  deformation  should  return  to  zero,  re- 
loading takes  place  along  the  original  loading  curve,  i.e.,  from 
1 to  2,  etc.  However,  if  the  deformation  increases  again  prior 
to  reaching  zero  the  member  will  reload  along  a straight  line 
from  the  point  where  reloading  starts  to  the  point  from  which 
unloc ding  had  started. 

3.2.2  Solution  of  Equations  of  Motion 

The  system  of  28  equations  is  solved  for  the  generalized 
accelerations  by  first  combining  the  vectors  on  the  right-hand 
side: 


[A]  {q}  = (B 1 } (91) 

where  (B'}  = {B}  + {P}  + {R}  + {Q} 
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Figure  30.  Example  of  Webbing  Force-Elongation  Curve. 


and  using  Gaussian  elimination  with  positioning  for  size: 


{q}  = [A]"1  {B} 


The  resulting  set  of  28  second-order  differential  equations  have 
the  general  form 

‘jj  = f j r ^2  ' “ ” ' ^28  ' ^3  ' ^2  ' * * * ' ^28^ 

q j (t=0)  = qAn  q4  (t=0)  = qAn  j = 1,  2,  ...,  28  (93 


These  equations  can  be  rewritten  as  56  first-order  equations 
having  the  general  form 


Yj  = fj  (t,  y1#  y2, 


• y28'  ql'  q2' 


'W 


q.  = y . 
j *3 


y-j  <t=0>  - 9j0  (t-0)  = qj0 


(94) 


Numerical  integration  of  this  set  of  equations  is  accomplished, 
using  the  Adams-Mculton  predictor-corrector  method  with  a vari- 
able step  size.  This  method  uses  the  difference  equations 

‘ yj,n  + I?  «***!.„  - + 37fj,n-2  '9f  J ,n-3)  (95) 
as  the  predictor  and 


(c) 

yj,n+l  * y 


24 


(9ffpl, . + 19f . „ - 5f . . + f . n ,) 

j , nt 1 j 3,nl  j , n 2 


(96) 


as  the  corrector.  Starting  values  are  provided  by  the  classical 
fourth-order  Runge-Kutta  method.  Input  data  includes  upper  and 
lower  error  bounds  for  the  solution.  Error  bounds  for  each  vari- 
able are  calculated  and  compared  at  each  step  with  the  difference 
between  the  predicted  value  yjp^  and  the  corrected  value  y.jc^  • If 
this  difference  exceeds  the  upper  bound  for  any  j , the  step  size 
is  halved.  If  this  difference  is  less  than  the  minimum  error 
bound  for  all  j and  for  three  successive  steps,  the  step  size  is 
doubled.  Halving  the  step  size  is  accomplished  by  interpolation 
of  past  data,  whereas,  doubling  is  effected  by  alternate  selection 
of  past  data.  The  solution  can  be  run  with  a fixed  step  size  by 
making  the  upper  and  lower  error  bounds  prohibitively  large  and 
small,  respectively,  or  by  using  equal  values  for  the  maximum 
and  minimum  step  size,  which  are  also  included  among  input  data. 
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3 . 3 PROGRAM  OUTPUT 


Output  data  consist  of  eight  blocks  of  information  that  are 
selected  for  printing  by  user  input.  The  data  include  time  his- 
tories of  the  variables,  which  are  simply  stored  during  the  solu- 
tion at  pre-determined  print  intervals  as  follows: 

1.  Occupant  segment  positions 

2.  Occupant  segment  velocities 

3.  Occupant  segment  accelerations 

4.  Restraint  system  loads  (tensile  loads  in  webbing  and 
resultant  normal  loads  on  pelvis  and  chest) 

5.  Seat  deflections  at  critical  points 

6.  Floor  reactions 

and  additional  information  as  follows: 

7.  Details  of  contact  between  the  occupant  and  the  air- 
craft interior 

8.  Injury  criteria 

The  last  two  blocks  of  output  data  will  be  discussed  in 
further  detail. 

3.3.1  Impact  Prediction 

For  prediction  of  impact  between  the  occupant  and  the  cock- 
pit interior,  23  surfaces  are  defined  on  the  body.  These  sur- 
faces were  illustrated  in  figure  15,  and  their  dimensions  dis- 
cussed in  section  2. 1.4. 4. 

The  distance  between  each  of  these  occupant  contact  surfaces 
and  the  aircraft  cockpit  surfaces  is  calculated,  during  execution 
of  the  program.  When  contact  occurs  between  an  occupant  surface 
and  a contact  surface,  the  time  and  relative  velocity  of  impact 
are  computed  and  stored  for  printing.  The  impact  conditions 
determined  in  this  way  can  be  used  in  evaluation  of  injury  poten- 
tial for  a given  cockpit  configuration. 
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3.3.2  Injury  Criteria 

The  injury  criteria  used  in  the  program  were  selected  as  the 
most  suitable  for  aircraft  crash  analysis.  The  dynamic  response 
index  (DRI)  is  used  to  determine  the  probability  of  spinal  injury 
due  to  a vertical  acceleration  parallel  to  the  spine  (reference 
10) . 

In  this  model,  the  response  of  the  body  to  acceleration 
parallel  to  the  spine  is  modeled  by  a single  lumped-mass,  damped- 
spring  system  as  shown  in  figure  31,  or,  in  other  words,  the  total 
body  mass  that  acts  on  the  vertebral  column  to  cause  deformation 
is  represented  by  the  single  mass.  In  general,  the  motion  of  the 
system  shown  in  figure  31  obeys  the  relationship 


2?oj 


dS 
n dt 


+ 


6 = Z 


(97) 


The  solution,  the  deflection  5,  is  representative  of  the  deforma- 
tion of  the  spine,  and  the  last  term  on  the  left-hand  side  of 
equation  (97) , divided  by  the  gravitational  acceleration,  is  the 
DRI.  The  properties  used  in  the  model  were  derived  from  tests  in- 
volving human  subjects  and  cadavers.  For  example,  the  spring 
stiffness  k was  determined  from  tests  of  human  cadaver  vertebral 
segments;  damping  ratios  were  determined  from  measurements  of 
mechanical  impedance  of  human  subjects  during  vibration  and  impact. 

In  the  occupant  model  used  in  program  SOM-LA,  it  is  assumed 
that  the  spring  in  figure  31  represents  the  lumbar  spine.  This 
is  a reasonable  assumption,  since  compression  fractures  that  occur 
in  vertical  impact  often  involve  this  part  of  the  spinal  column. 
Therefore,  the  mass  m in  figure  31  is  the  sum  of  segments  2,  3, 

4,  5,  6,  and  7.  The  input  acceleration  is  the  component  of  the 
acceleration  of  segment  1 in  the  z^-direction: 


*1  ' *1  T13  + *1  T23  + *1  T33 


(98) 
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mass  (lb-sec  /in.) 
deflection  (in.) 
damping  ratio 
stiffness  (lb/in.) 
acceleration  input 
(in. /sec2) 


*DRI 
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*Dynamic  Response  Index 


wn  = natural  frequency  of 

the  analog  = 
(rad/sec)  \ /m 

g = 386  in. /sec2 


Figure  31.  Model  Used  for  Prediction  of  Spinal 
Injury  (From  Reference  10) . 


so  that  the  equation  to  be  solved  is 


6 + + u>n26  = z 


(99) 


where  £ = 0.224 


to  =52.9  rad/sec. 
n 


The  DRI  is  calculated  at  each  step  by 


DRI  = to  S/g 
n 


U00) 


and  the  maximum  value  is  stored  for  output. 

Another  injury  criterion  is  the  Severity  Index  (SI)  (refer- 
ence 11) , which  is  calculated  for  the  head  and  chest 
according  to 


tf 

SI  = / a ' dt 

fco 
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(101) 


where  a = acceleration  as  a function  of  time 


n = weighting  factor  (n>l) 
t = time 


The  severity  index  has  been  validated  for  frontal  impacts  of  the 
head-face  with  n=2.5.  Although  such  a validation  has  not  been 
performed  for  other  parts  of  the  body,  the  SI  is  calculated  for 
the  chest  as  it  may  serve  a useful  function  in  determining 
relative  levels  of  injury  potential  due  to  an  acceleration : 
environment. 

Finally,  the  Head  Injury  Criterion  (HIC)  contained  in  Federal 
Motor  Vehicle  Safety  Standard  20 J is  calculated  according  zo 


(102) 


where  a is  the  resultant  head  acceleration  and  4 and  t2  are  any 
two  points  in  time  during  the  crash  event. 
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4 . 0 SUMMARY 


A three-dimensional  mathematical  model  of  an  aircraft  seat, 
occupant,  and  restraint  system  has  been  developed  as  an  aid  to 
the  development  of  crashworthy  seats  and  restraint  systems  for 
general  aviation  aircraft.  The  occupant  model  consists  of  eleven 
rigid  mass  segments  whose  dimensions  and  inertial  properties 
have  been  determined  from  studies  of  human  body  anthropometry 
and  kinematics.  The  seat  model  is  made  up  of  beam  and  membrane 
elements  with  provision  for  simulating  plastic  behavior  by  the 
introduction  of  plastic  hinges  in  the  beams. 

A user-oriented  computer  program  called  Seat  Occupant  Model- 
Light  Aircraft  (SOM-LA)  based  on  the  three-dimensional  model  has 
been  developed  for  use  by  engineers  concerned  with  design  and 
analysis  of  general  aviation  seats  and  restraint  systems.  De- 
tailed descriptions  of  both  are  used  as  input  data.  The  response 
of  the  seat  and  occupant,  restraint  system  loads,  and  various  in- 
jury criteria  are  predicted  for  any  given  set  of  crash  conditions. 

Results  of  the  computer  program  validation  will  be  covered 
in  a report  entitled  "Validation  of  the  SOM-LA  Computer  Program." 
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APPENDIX  A 

OCCUPANT  SEGMENT  POSITION  AND 
VELOCITY  COMPONENTS 


A-l 


A. 1 SEGMENT  POSITION  COMPONENTS 

Referring  to  figure  13,  which  has  been  repeated  here  as 
figure  A-l  for  convenience,  and  using  Pn  = Ln  ” Pn  t*16  absolute 
position  of  the  mass  center  of  each  body  segment  is  given  below. 
The  elements  of  the  transformation  matrices  [T11]  are  functions 
of  the  generalized  coordinates,  as  given- by  equations  (6),  (9), 
and  (10)  . 


Segment  1: 


(X^,  , Z^) , the  coordinates  of  the  reference  point  on  tne 


Segment  2 : 

X. 


Segment  3: 
X.  ‘ 


0 

0 

P. 


<qr 

q3> 

- -> 

• 

0 | 

T2 

0 

- 

p2  J 

- - 

0 ] 

T2 

0 

L2) 

T3 

Lc 

0 

°3 

Segment  4: 


A-2 


+ 


Segment  10: 


Segment  11: 


| 


A- 5 


I 

i 


A. 2 SEGMENT  VELOCITIES 


The  velocity  components  for  the  mass  center  of  each  segment 
are  given  below.  Here,  the  symmetry  of  the  body  has  been  used 
for  convenience  in  combining  similar  terms  at  a later  point.  By 
symmetry , 
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Z9  = Zi  - PxT33  - LflT32  + LgT31  - Pg^33 


• • • 1 *1  *10 
X10  " X1  * P1T13  + LHT12  + P8T11 


*10  “ X1  ~ Pl*23  + V22  + p8*21 


• • • 1 • 1 • IQ 

Z10  Z1  " P1T33  + LHT32  + P8T31 


• • • 1 • 1 • IQ  • 11 

X11  = X1  “ P1T13  + LKT12  + L8T11  “ P9T13 


*11  X1  " pl*23  + LH^22  + L8*21  ' p9*23 


• • *1  *1  • 10  • 11 
Z11  = Z1  P1T33  + LHT32  + L8T31  “ P9T33 


A- 8 


APPENDIX  B 

OCCUPANT  KINETIC  ENERGY 
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Occupant  Kinetic  Energy 

The  system  kinetic  energy  T is  written  below  as  a function 
of  the  generalized  coordinates.  Prior  to  the  kinetic  energy  a 
set  of  50  constants,  functions  of  body  properties,  are  defined. 
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z10  = Z1  P1T33  + LHT32  + P8T31 


10 


Z11  = Z1  “ P1T331  + V321  + L8T31  ° " P9T33 


C-.1 


i 


APPENDIX  D 
[A],  (B) 


« 

! 

I 


D-l 


■ ■ zsmsmm^mrnBgmm 


Ef 


Below  are  listed  the  non-zero  elements  of  the  matrix  [A] 
and  the  vector  (B>,  derived  from  the  kinetic  energy  derivatives 
Because  [A]  is  symmetric,  only  the  upper  triangular  part  is 
given  to  save  space,  i.e. 


I 


A^j  where  j _<  i 


A (1,1) 


= C, 


I 


A (1,4)  = CuT23j 


A (1,5)  = CnT11  cos  <J>1 


A (1,6)  = - CUT12J 


A (1,7)  = - C12T23 


A (1,8)  - C12T1]L  col  <}>2 


A (1,9)  = - C12T12' 


a (1.10)  = - c13t233  - c14t213 


A (1,11)  = (C13T333  + C14T133)  cos  *. 


A (1,12)  = - C1;jT12' 


A (1.33)  = C15T234  - C16T215 


A (1,14)  = (C16T315  - C15T334)  cos  ^ 


A (1,15)  = (C15  + cos  a5)  T12 


A (1,16)  = C16T13‘ 


D-2 


A (1,17)  - C15T236  - C16T21 

A (1,18)  = (C16T317  - C15T336)  cos  »6 

A (1,19)  = (C15  + C16  cos  a?)  T126 

A (1,20)  = C16T137 

A (1,21)  = - C17T218  + C18T239 

A (1,22)  = - (C17  sin  6g  + C18T339)  cor 

8 

A (1,23)  = C18T12  sin  a& 

A (1,24)  - - C18Tn9 
A (1,25)  = - C17T2110  + C18T23U 
A (1,26)  = - (C17  sin  S10  + Cj^jTjj11)  cos  *10 
A (1,27)  = C18T12  sin  a31 
A (1,28)  = - C18Tnn 
A (2,2)  =■  C3 

A (2,4)  = C^j1 
A (2,5)  = Cj^T^1  cos  4^ 

A (2,6)  = - Ci;lT22 
A (2,7)  = C12T132 


D-3 


A (2,8)  = c12T21  COS  ^2 


A (2,9) 


" C12T22 


A (2,10)  = 


A (2,11) 


13T133  + C14T113 


(C13T333  + C14T313)  Si”  *3 


A (2,12)  = 


" C13T22 


A (2,13)  = 


A (2,14)  = 


- C15T134  + C16T115 
<C16T315  ' C15T334)  Sin  *4 


A (2,15)  = (C15  + Clg  cos  a5)  T22 


A (2,16)  = 


C16T23' 


A (2,17)  = - Cl5Tu6  + C16Tu7 
A (2.18)  = (C16T317  - C15T336)  sin  *6 


A (2,19)  = (C15  + Clg  cos  a?)  T22 


A (2,20) 


- C T 
^16  23 


A (2,21)  - Ci?T118  - C18T13 


A (2,22)  = - (C17  sin  0g  + C]_8T33  ) sin  ^8 


A (2,23) 


= C18T228  Sin  “9 


A (2,24) 


= “ C18T21 


D-4 


y •*  >1  ■*  y 


Kj?'_ 

TS 

I 


I 

I 


I 


r. 


K t 


l 


C T ^ - C T ^ 
7^11  1 


'18  13 


11, 


A (2,25)  = 

A (2,26)  = - (C1?  sin  eiQ  + C18T33^)  sin  H>1Q 
A (2,27)  = 


C18T2210  Sin  all 


~ Cj.8T21 


11 


= C. 


A (2,28) 

A (3,3) 

A (3,5) 

A (3,6) 

A (3,8) 

A (3,9) 

A (3,11) 

A (3,12) 

A (3,14)  = C16(sin  04  cos  <}>4  cos  a5  - cos  04 


C11T31  COS  *1 


_ r T 
'“ll  32 


C12T31  COS  *2 


- C T 
''12i32 


C13T31  cos  *3 


-C  T 
13*32 


C. . cos  0O 
14  3 


- C15T31 


(Clc  + Clc  cos  a5)  T32 


15  *16 

5 


r T 
'16  33 


A (3,15) 

A (3,16) 

A (3,18)  = Clg  (sin  06  cos  <t>g  cos  a?  - cos  0( 

_ r T 
L15  31 


D-5 


sin  a,.) 
cos  <f>4 

. sin  a^/ 
’ cos  <j>g 


A (3,19)  = (C15  + Clg  cos  a?)  T32 


A (3,20)  - C16T33‘ 


A (3,22)  = - C1?  cos  0 


- C,  Q (cos  0 q cos  aQ  - sin  0Q  cos  sin  aQ) 


A (3,23)  = C18T328  sin  ag 


A (3,24)  = - C18T31' 


A (3,26)  = - C1?  cos  610 


- C18  (cos  01Q  cos  - sin  01Q  cos  <j>10  sin  a^) 


A (3,27)  = C18T3210  sin  <*u 


A (3,28)  = - ci8T31 

A (4,4)  = (C2  + Cg)  sin2  + C37  sin2  0.^ 

2 2 2 
+ (C2  sin  0^^  + Cg  + C37  cos  0^  cos  <j>^ 


A (4,3)  - (C2  ~ Cg  - C37  + C38)  T32  cos  ({^ 


A (4,6)  = (C2  + Cg  + C3?)  T31 

A (4,7)  = C19  (T13XT132  + T^1*^2) 

A (4,8)  = C19  cos  <f>2  (T13  T21  - T23  T1a  ) 


A (4,9)  ci9(T2i  T12  " " T13  T22  } 


* »%**■  wS^W'S'TcT^ 


'.  * r , 
A V.  '' 


A (4,10)  - C2|)  (T^T^3  + Tjj'-Tjj3) 


1„  3 


1„  3, 


+ C21  (T13  T11  + T23  T2L  * 

’ T * 

*20A33 

i m 

'20  33 


A (4,11)  = (ConT_3  + C21T313)  T13X  sin  ip3 
" <C,«T„3  + C21T313)  T231  COS  ^ 


A (4,12)  = C2Q  lt23\23  ~ T131T223) 

A (4,13)  = - C25  (T231T234  + T13XT134) 

+ C28  ( T23lT215  + T131t115) 

A (4,14)  = (C25T33  - C2gT31  ) (T23  cos  i|>4  - T13 


A (4,15)  = (C23  + C2g  cos  a,.)  (T13  T22  - T23  T12 

A (4,16)  = C2g  (T13  T23  - T23  T13  ) 

A (4,17)  = - C25  (T13  + T23  T23  ) 


+ C28  (T23lT217  + T13±T117) 

A (4,18)  = <C25T336  - C28T317)  (T^1  cos  t|>fi  - T13X 


23  ww“  T6  13 

lm  6 „ 1„  6 


A (4,19)  - (C25  + C2g  cos  a?)  (T13  T22  - T23  T12 


28 

1 7 

(T13  T2  3 

" T23 

34 

(1  XT  9 
v 13  13 

+ T2  3 

35 

1 9 

(T  "T  ^ 
' 1 12  1 13 

+ T22 

XT  7) 
13  1 

XT  9) 
X2 3 J 

1„  9. 


sin  ij>4) 

) 


sin  <^g) 

) 


D-7 


Lr  <" /*"  • <*-]  * - „KvV».  * • > - V " rv7  ■ f v>' *• 


I 


K 


1 


i 


i 


_ C32  (T13  T11 


•IsSl8  + T23lT-8) 


r at  a) 
r23  T21  1 

- S3  + T22lT218) 

A (4,22)  - (C34T339  - C32T318)  (T13Ssin  t8  - COS  V 

+ (C3'T339  " C33T318)  (Tx21  sin  ^8  " T',° 

A (4,23) 

*>]  si 


[C34  (T23  T12 
+ C 


'221  cos  V 

, 1„  8 m 1»  8 

•23  T12  T13 


• XT  °) 
13  *22  ' 

18-1. 

:35  ‘T22  T12 


;121t228)J  Sin 


9 

A (4,24)  = C34  (T131t2i9  " T23iTH  * 

+ C35  (T12  T21 


. 1*  9 _ t XT  9) 

12  T21  T22  11 


A (4,25) 




1 11  ml, 

=34  (T'3  T13  + T23 


i3U  + T23lT23U> 
=35  <T12lT13U  + T221' 


r22lT23U 


) 


c (T  V,10  + T,,1^,10 
C32  ' 1 13  11  23  21 


1 10  V 

r23  T21 

+ =33  (T121t1110  + T221t21W) 

(4.26)  = (C34T33U  - =32T3110)  si"  *10  " ^ °°S  ♦l0) 

- (=35T33U-  =33T3110)  sin  *10  ‘ ^ 

(4.27)  = [c35  (T121t22J°  " T221t1210) 

c (T  S 10 
C34  U13  i22 


^lO 

r221  C0S  V 


V°  - -231-121°>] 

,28)  = C34  (T131T2111  ’ T231tU:I'1) 


sin  a 


11 


D-8 


1 

- C35 

(T121t21U  - T22lTll 

1 

1 

A 

(5,5)  = (C2 

+ C3g)  cos2  + (Cg 

! 

t 

% 

A 

(5,7)  = Cw 

■\  7 12 

(T21~T13  “ T11  T23 

! 

• 

A 

(5,8)  = Cig 

cos  ^ cos  4>2  (t111t 

m 


m 


Kv  : 

K,  ' 


:6  + c37)  sin2 


2 12 
i z 4-  T XT  £ 
11  a21  i21 


+ T311t312) 

A (5,9)  = - C19  cos  4»1  (Ti;lXT122 

A (5,10) 


2 12  12 
1 c + T aT  4 + T XT  ) 
12  21  22  31  32  ; 


[<C20T133  + C21T113)  T213 
(C20T2  3 


A (5,11) 


:20T23  + C21T21  J 'll  ] COS  *1 

C. 


[C20  COS  ^3  (T111t113  + T21  T21  + T31 

'21  {(T11 


i XT  3) 
31  A31  1 


r211  sin  ij>3)  sin 

+ t312  oos  63j] 


'll1  oos  *3  + T21l 

9 3) J CCS 


0. 


A (5,12) 
A (5,13) 


“ C20  COS  ^1  (T11  T12  + T21  T22  + T31  T32  J 

[C25  {T11  T23  " T21  T13  } 


2 (T  3 

28  U11  A21 


T21lrll5)J  COS  +1 

'25  (TU  T11  + T21  T21  + T31  T31  1 cos 


Coo  T21 

A (5,14)  = - C, 

< 

[(Tu1  oos  + t21X  sin  *4)  T315 

3.  COS  <j).  cos 
4 4 


+ C, 


'2  8 


+ T-^^sin  6, 


a 


:c  - cos  0 . sin  a^)  cos 
5 4 d J 


D-9 


A (5/15)  = 
+ 

A (5/16)  = 
A (5/17)  = 

A (5.18)  = 
+ 

A (5,19)  - 

-} 

A (5,20)  = 
A (5,21)  = 

A (5,22) 


{C25  + °28 
T311t324)  COS  *1 


cos  ct,.)  (T^^  T^2 


T12  21  22 


COD  cos  <tx 


'28 

[C25  (T11  T2 3 


15  1 ! 

. T13  + T21  T23 


*1  (T11  T13 
Im  ^ _ <T 

. T2  3 21  A13  * 


IfTi  5 , m lrp 

. T23  31  X33 


') 


C28  (TlllT217  - T21lTU7)]  °Ci  *1 

1 cos  + T211  sin  \P6)  T31 

<p  cos  a<7 


L11 


“28 

+ I C28  [(T3 

T311  (sin  06  cos  <p6 

1T316)  cos  <D6J 


- cos  6-  sin  a_ 
6 ! 


a J 


1«P  6 4.  7 ^ 

1 T21  h 31  1 31 


C25  (T111t116  + T21lT216 

1_  6 


cos 


1T326)  COS  4) 


; (C25  + C28 
• T31lT32 
:28 
(C32 
(C34 


C2 8 cos  a7)  (T111t126  1 T22  \-.2 


= C28  cos  (t>1  (T111T13'  + T211t237  + T31  T33  } 

lm  ® _ <p  ^ 

T21  ^21  Ali 


cos  <^>1  + C33  sin  4>-|_)  t2i 


cos  4^  + C35 


= (C32  cos  4>1  + c 


33 


sin  *x)  (Tj^^Tjs9  - T211T1 
sin  <j>1)  [(Tu 


+ T21X  sin  ^8)  sin  eg  + T31' 


1 COS  0g 


1 COS  l|/g 

.] 


+ (c34 


cos  <t>1  + C35 


+ ^l1  sin  V-  T33 
- sin  Og  cos  4> g s^n  ag 


sin  <j>^)  cos  08 

9 

+ T31 


) 


1 ( cos  0g  COS  ag 


D-10 


A (5,23) 


A (5,24) 


A (5,25) 


A (5,26) 


A (5,27) 


A (5,28) 


- <C34  COS  4^  + c35  sin  *!> 


T21lT228  + T31lT328)  sin  “< 


(C34  cos  ^ + C35  sin  41)  (Tjj1^9 


19  19 

T21  T21  + T31  T31  ) 


(C32 

cos 

*i  - 

C33 

sin 

♦i>  (tii1t2110 

- T 1T  10 
zl  L11 

(C34 

cos 

*i  - 

C35 

sin 

♦l>  <T111t23U 

- T XT  11 
l21  A13 

(C32 

cos 

*i  - 

C33 

sin 

61)  [(THl  COS 

910 

T 1 
i21 

sin 

*10> 

sin 

910 

+ T311  cos  0 10 

] 

(C34 

cos 

♦l  * 

C35 

sin 

♦i>  [<Tii1-00S 

*10 

T X 
i21 

sin 

*10  > 

T33 

11  + 

T311  (cos  e1(J  ■ 

cos 

sin 

610  ' 

COS  <f> 

io  sin  a 

n>] 

(C35 

sin 

♦l  " 

C34 

cos 

*i>  (tii1ti210 

► T211t2210  + T311t3210>  sin  “ll 


lm  11 


(C34  cos  <j>1  - C35  sin  <|>1)  (Ti;l  Tu 


. m l|p  11  l m ^*33 

f X21  i21  31  “31  1 


D-ll 


/ 


¥ 

F 

I 

£~ 

& 

% 

I 


jt 

ir 


| 

$ 

i 

a 

t 


f 

I 

I 


jr 

l 


i 


A (6,6)  = C2  + Cg  + C3? 

A (6,7)  = C19  (T121T233  - T221T132) 

A (6,8)  = - C19  COS  <J>2  (Ti2  T.^  + T22  T21  + T32  T31  ) 

A (6,9)  = C19  (T12  T12  + T22  T22  + T32  T32  ) 

A (6,10)  = C2()  (T121T233  - T221T133) 

+ C21  (T12  T21  " T22  T11  1 

A (6,11)  = - C2Q  cos  4> 3 (T12  Tx1  + T22  T21  + T32  T31  ) 

+ C21  [^T121  COS  ^3  + T221  sin  ^3^  sin  03 

+ cos  e3] 

A (6,12)  = C2Q  (Tu1T123  + T221T223  + 

A (6,13)  = C25  (T134  T223  - T234  T123) 

C28  <T11  T22  " T21  T12  * 

A (6,14)  = C25  cos  cf4  (T12  T1;l  + T22  T21  + Tj2  T31  ) 

' C2  8 LT315  <T121  cos  *4  + sin  V 

+ T^1  (sin  0^  cos  <j>4  cos  a5  - cos  04  sin  a5)J 

A (6,15)  = - (C25  + C28  cos  a&)  (T121T124 

i 4 14 

+ T T + T T ) 

22  22  32  32  ' 


D-12 


A (6,16)  = - C28  (Tu1T135  + + T321T335) 


6m  1 


6m  1, 


A (6 , 17)  = C25  (T13  T22  - T,,  T19  ) 


23  12 


C28  *T11  T22  “ T21  T12  J 

A (6 ,18)  = C25  cos  <j»g  {TyV  + T221T216  + T321T316) 

" C28  [T31?  (T121  COS  ^6  + T22X  sin  V 
+ T.^1  (sin  0g  cos  4>g  cos  a?  - cos  0g  sin  a7)j 

A (6,19)  = - (C25  4 C2g  cos  a7)  (T121T126 
+ T221t22  + T32  T32  ) 

A (6,20)  = - C2g  (T12  4 T22  T23  4 T32  T33  ) 


A (6,21)  - - C32  (T12  T21  - T22  T13  ) 


+ C34  (T12  T23  " T22  T13  ) 


+ C33  (T13  T21  ' T2 3 T11  ) 
C35  (T13  T2  3 " T2 3 T13  ) 


A (6,22)  = - (C32  sin  0g  4 C^T^9)  (T.^1  cos  i|»g 


+ T221  sin  iP8) 


D-13 


A (6,23) 


A (6,24) 


A (6,25) 


A (6,26) 


[c32  cos  0g  + C34  (cos  0g  cos  ag 

sin  0g  cos  4»g  sin  ag)j  T.^1  + (C33  sin  0( 


' C35T339)  ‘’’IS1  008  *8  + sin  *8’ 


[c  , cos  0O  + C,c  (cos  0O  cos  a 


8 35 


8 


sin  0g  cos  4>8  sin  ag)]  T.^1 


[C34  (T121t128  + T221t22  + T32  T32  } 

C35  {T13  rj?12  + T23  T22  + T33  T32  0 


" C34  (T121t119  + T22  T21  + T32  T31  ) 


f C35  (T131t119  + T23  T21  + T33  T31  * 


- C (T  - T ^0 

C32  U12  i21  l22 


+ c34  (T12lT23n  * T221t13U) 


_ c (T  - T ^) 

C33  T13  A21  1 23  ill  1 


* C35  (T13lT23n  - T231t13U) 


- ' C32  COS  *10  + T221  Sl"  *10>  Sin  9 


10 


+ T32  cos  01O 


] 


D-14 


• ■3t3^Sfg^^SiXairt--r 


11 


“ C34  [(T12  COS  ^10  + T22  sin  ^10)  T33 


+ ,»,321  (cos  610.  cos  ai;L  - sin  01Q  cos  <J>1()  sin  ai;L)J 


- C33  [(T^1  cos  *1Q  + T^1  sin  *10  ) sin  9 
] - c35  [^n1  cos  ho 


10 


+ T33  cos  e10j  - c35 


+ T^1  sir.  *10>  T,/1  TJ31  ( cos  «10  cos  «u 


10'  33 


33 


10 


- sin  01O  cos  4>10  sin  a^)! 


A (6,27)  = [c34  llj/ljj10  + T221t2210  + T321t3210) 

+ C35  (T131t1210  + T23lT2210  + T331t3210>]  sin  °11 
A (6,28)  = - CJ4  (T121Tu11  + T./t^11  6 T321T3111) 


- C35  (T131tU11  + T23lT21U  + 


2 2 2 

A (7,7)  = C3  (sin  02  cos  <t>2  + sin  <j>2) 


2 2 2 
+ C5  (sin  02  sin  <j>2  + cos  <J>2) 


+ C39  sin2  02  + C4Q  cos2  $2 


A (7,8)  = (C3  - C5)  T32  cos  <t>2 


•utessa 


D-15 


A (7,9)  = (C3  + C5  + C39)  T31" 


A (7,10)  = C22  (T13  T13  + T23  T23  ) 


+ C23  (T13  T11  + T23  T21  } 


A (7,11)  = (C22T333  + C23T313)  (t132  sin  t|>3  - T233  cos  4^) 


A (7,12)  C22  (T12  T23  - T22  T13  ) 


A (7,13)  = - C2g  (T132Tn4  + T232T234) 


+ C27  {T122t134  + T222t234) 


+ c29  (T132T11  + t23  t21  ) 


- c30  {T122t115  + T222t215) 

A (7,14)  = (C26T334  - C2gT315)  (T232  cos  *4  - T132  sin  *4> 

- (C27T339  ~ C30T315)  (T222  COS  *4  ' T122  sin  V 


A (7,15)  - (C26  + C2g  cos  a5)  (T13  T22  - T23  T12  ) 


V '13  22 


2„  4 


2„  4. 


" (C27  + C30  COS  a5)  (T12  T22  “ T22  T12  ) 


A (7,16)  - C2g  (T13  T23  - T23  T13") 


" C30  (T12  T23  ” T22  T13  ] 


D-16 


It 


I 


K 


2„  3 


2„  3 . m 2„  3. 


A (8,12)  = - C22  cos  <(>2  (T^  T12  + T21  T22  + T31  T32  ) 

2 4 2 4 

A (8,13)  = (C2g  cos  4>2  - C 27  sin  <t>2)  (T.^  T23  - T21  T13  ) 


- (C 29  cos  <(>0  - C^n  sin  4>0)  ( 


2 30 


T 2T„,5  - T,,2!.5) 


r2'  '"11  "21  21  ‘11 


A (8,14)  = - (C26  cos  <p2  - C27  sin  *2)  (T112T114  + T^T^4 


+ T312t314)  COS  ^4 


+ (c29  cos  (J>2  - C3Q  sin  <J>2)  [(Tn2  cos  i|>4 


+ T,,2  sin  ^4)  T315  + T312  (sin  04  cos  <j>4  cos  a5 


- cos  04  sin  a^)J 


A (8, IS)  = [(C2g  + C29  cos  a5)  cos  <j>2  - (C27 

+ C3Q  cos  a5)  sin  4>2]  (Tn2T124  + T212T 


2™,  4 

22 


+ t3i2t324) 


A (8,16)  = (C29  cos  (j>2  - C3U  sin  <J>2)  (Tn2T135 
+ T21  T23  + T31  T33  J 
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2 6 2 6 

A ( 8,17)  = (C26  cos  *2  + C2?  sin  <f>2)  (Tu  T23  - Tn  T±3  ) 

- (C2g  cos4»2  + C3Q  sin  <J>2)  (T^  T21  - T21  ) 


A (8,18)  = -(Coc  cos  <P2  + C27  sin  <}>,)  (Tn  T1T  + T21  T 


2 x 11  11  21  21 


+ T312T316)  cos  <j>g  + (C29  cos  <{»2 

+ C30  sin  <02)  [(Tn2  cos  + T212  sin  T^7 

+ T312  ( sin  6g  cos  <p6  cos  a?  - cos  0g  sin  a7>J 

A (8,19)  = [<C26  + C29  cos  a?)  cos  <j>2  +(C27 

+ C3Q  cos  a?)  sin  <|>2J  (Tu  T12  + T21  T22 


+ t31  t32  ) 


2„  7 


2m  7 


A (8,20)  = (C2g  cos  <J>2  + C3Q  sin  <J>2)  (Tn  T13  + T21  T23 


+ T 2T  7) 

+ in  i33  ; 


A (9,9)  = C3  + C5  + C39 


A (9,10)  - C22  (T12  T23  - T22  T13  ) 


f C23  (T12  T21  " T22  T11  * 
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A (9,15) 


+ 


A (9,16) 


A (9,17) 


A (9,18) 


' (C26  + C29  COS  “51  (T122t124  + T222t224 


- T 2T  ^) 
a32  X32  } 


(C27  + C30  “s  “s’  (T132t124  + T232t224 


H T 2T  ^ ) 
H T33  i32  ) 


- C29  <T122t135  + T222t235  + T322t335> 


- C30(T132t135  + T232t235  + T332t335> 


" ' c26  (T122t23  " T22  T13  ’ 


4-  r (T  ^ - T ^ T ^) 

+ C27  T13  T 6 23  i13  ' 


+ C29  (T122t217  ' T222t117’ 


2 7 „ 2_  7. 

" C30  *13  T21  " T23  T11  ’ 


- [C26  <*12  V + *222*216  + *32‘*3l’> 

' *"27  (T122t116  + T222*216  + T322t316’]  “°S  *« 


2„  6, 


- C 


2 cos  ^ + T002  sin  ’l»6)  Tn 


29  [(T12  r 6 • x22 

+ T322  (sin  Qg  cos  <t>6  cos  cc  7 - cos  9g  sin  a7)j 

+ C30  [(T132  COS  *6  + T222  Sin  V T317 


+ t33  (sin  06  cos  <|>g  cos  a. 


- cos  0,  sin  3tJ 
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A (9,19)  = - (C26  + C29  cos  a?)  (T122T126  + T 222T22 6 
+ T322t326)  + (C27  + C30  COS  a7)  {T132t126 


+ T23  T22  + T33  T33  * 


A (9,20)  = - C2g  (T12  T13  + T22  T23  + T32  T33  ) 


+ C30  (T132t137  + T232t23?  + 


2 2 .2 

A (10,10)  = C13P3  (sin  03  cos  <j>3  + sin  <fr3) 


7 2 2 

+ c41  (sin  03  + cos  03  sin  <t>3) 


2 2 2 
+ C42  cos  9^  cos  <t>3  + C4  cos  9. 


+ 2 C24T333  sin  93 

A (10,11)  = [(C13P3  + C41  - C42)  T333  + C24T313]  sin  ^ 


A (10,12)  = - (C13P3  + C41)  sin  03  - C24T33' 


(C13D?  + C41)  T31  " C24 T 3 3 


A (11,11)  = (C13P3  + C41)  cos2  4>3  + C42  sin2  d>3  + C4 


A (11,12)  = C24  sin  4>3 


A (12,12)  = C13P3  + C41 


D-22 


'31 


t5)  (sin2  04 


[< 


2 

*4  ♦ 

sin' 

cos 

COS  <t»4 

COS 

a5,: 

4 2 

4 

1 

) + 

T32 

> J 

+ C44  lT334>2  + C46  <T335>2  + C45  O 


+ 2 C31T31  T33  sln  “5  + C46  <T32  * 


A (13,14) 


-k 


+ C.-  - C..  + 2 C,.  cos  ac 
43  44  31  5 


+ (Cg  - C45  + C4g)  cos 


a5j  cos  04  sin  <j>4  cos  <j>4 


+ j^(Cg  - C45  + C4g)  cos  a. 


+ sin  04  sin  <J>4  sin  a,. 


. 2 2 

- Cg  sin  \j>4  cos  ij>4  sm  04  cos  04  sin  <J>4  cos  a5 


A (13,15)  = (C_  + C31  cos  a5  + C43)  Tgi 


+ (C8  + C46)  T33  COS  a5  + C31T33 


+ C45T31  Sin  ?5 


A (13,16)  = (C46  - Cg  - C31  cos  ag)  T 


32 
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A (14,14) 


A (14,15) 
A (14,16) 
A (15,15) 


A (16,16) 
A (17,17) 


- (C?  + C43  + C4g  + Cg  cos  a5 


+ 2 C3^  cos  a,.)  cos 


+ (C44  + C45  cos2  a5  + c46  sin2  a5>  sin2  *4 


+ cg  sin2  a5 


~ “ (Cg  - C45  + C4g)  cos  a7  + C31  | sin  cbg  sin  a 


(C46  ' C8  ' C31  COS  a5>  COS  ^4 


2 

C_  + C..,  + 2 C-,  cos  ac  + C . c sin  ac 
7 43  31  5 45  5 


+ (Cg  ^4g)  COS  Otj. 


* C8  + C46 


[2  2 
(C7  + 2 C31  cos  a?  + Cg  cos  a7)  cos  <f> ^ 

2 2 1 . 2 
+ C43  + C45  sin  a7  + c46  cos  a7J  sin  6g 


+ (C44  + C46  cos2  a,  + C46  sin2  a?)  (T336) 


+ (C«  + C46>  (T32  > 


2 2 
+ (C7  + Cg  cos  a7  + 2 C31  cos  a 7)  sm  <j>  g 
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[< 


2 I (Cg  C45  + C4g)  cos  *7 


16  2 2 
+ C.ji  sin  sin  ct^  + Cg  cos  6g  sin  <*7 


A (17 , 


[c7  + 

C43  " 

C44  + 2 C31  COS  a7 

(C8  - 

C.5  + 

c46)  cos2  a7 j T326 

COS  <J>g 

i 1 

n 

00 

- c45  ^ 

’ C4g)  cos  a? 

C3l] 

sin  0- 
0 

sin  <j>g  sin  a? 

Cg  sin  cos  sin  0g  cos  0g 

2 2 
sin  <j)g  cos  a7 

(C?  + 

C43  + 

C3i  308  a7^  T3i 

1 1 

n 

00 

+ C46) 

cos  a?  + C31]  T337 

+ C45T31  sin  a7 

(C46 

- C8  - 

C31  cos  a?)  T32 

(C?  + 

C43  + 

2 C31  cos  a7  + Cg 

2 

cos  a7 

+ c46)  cos2  + (C4 4 + c45  cos2  <*7 


2 2 2 
+ C.,  sin  a_)  sin  <(>..  + C_  sin  «. 

4o  / DO  / 

A (18,19)  = - j^(Cg  - C45  + C46)  cos  a7  + C31j  sin  sin  a7 
A (18,20)  = (C,,-  - CQ  - C_.  cos  a_)  cos  dr. 

40  O J I / t 
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A (19,19)  = C?  + C43  + 2 C31  cos  a? 

+ (C8  + C46)  COs2  “7  + C45  sin2  a7 
A (20,20)  = Cg  + C46 

A (21,21)  = [C9  + C48  + (C49  cos2  + C10)  cos2  ag 
+ C49  sin2  <t>8  + 2 C36  cos  ctg 
+ C50  cos2  4>g  sin2  a9J  cos2  0g 
+ [c47  + (C1Q  cos2  0>8  + C49)  sin2  a9 
+ C5Q  cos2  ag]  sin2  0g  + C1Q  sin2  <j>8  sin2  ag 

- 2 [(cio  _ C49  " C50}  COS  a9 
+ c36]  sin  08  T338  sin  a9 

A (21,22)  = (C1Q  + C49  - C5Q)  T339 

+ C36  sin  0 8 sin  <j>8  sin  ag 

A (21,23)  = [<C10  + C49  - C5Q)  cos  a9  + C3g]  T33®  sin  ag 

' [(C10  + C49}  Sin2  a9  + C50  COs2  a9  + C4 7]  sin  08 

A (21,24)  = (C10  + C49  + C36  cos  a9)  T328 
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A (22,22)  = Cg  C4g  + 2 C3Q  cos  ag 

+ (C10  + C49  sin2  V COs2  a9  + C49  cos2  4*8 
+ (C1Q  cos2  cj>g  + C50  sin2  <08)  sin2  ag 

A (22,23)  = - [(C10  + C49  - C5Q)  cos  ag  + C3gJ  sin  4>g  sin  ag 

A (22,24)  = (C1Q  + C49  + C36  cos  ag)  cos  <J>g 

A (23,23)  = (C^8  + C4g)sin  ag  + Cgg  cos  ag  + C 47 

A (24,24)  = C1Q  + C4g 

A (25,25)  = |^Cg  + C48  + (C10  + C4g  cos2  <!>10)  cos2  an 
+ 2 c36  cos  au  + C49  sin2  +J0 
+ C5Q  cos2  ♦10  sin2  ,u]  cos2  e10 

+ [(C10  cos2  +10  + C49>  sin2  “ll  + C4? 

+ C50  COs2  “ll]  sin2  °10  + C10  S“2  +10  sin2  “11 
+ 2 [<C10  + C49  " C5o’  °°s  “ll 
+ C36 j T3110t3310  sin  “ll 
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f 

| 


A (25 ,26)  = [(C10  - C49  - C5Q)  T331 

+ C36  Sin  \o]  Sin  ^10  Sin  “ll 
A (25,27)  = £(C,0  + C49  - Cr>0)  cos  + C3gj  T33  sin 

+ [{C10  + C49J  sin  all  + C50  °°S  all  + C4 7]  T31 
A (25,28)  = (C1Q  + C49  + C3g  cos  a^)  m32 

A (26,26)  = C9  + C4g  + 2 C3g  cos 

+ (C1Q  cos2  4>10  + C5Q  sin2  4>10)  sin2  <*u 

+ (C10  + C49  Sin2  'V  COs2  “ll  + C49  COs2  ^10 
A (26,27)  - - [(C1Q  + C49  - C5Q)  cos  «u  + C3g]  sin  <1>10  sin  au 

A (26,28)  = (C1Q  + C49  + C3g  cos  au)  cos  4>10 

A (27,27)  = (C1Q  + C49)  sin2  an  + C5Q  cos2  au  + C 4? 

A (28,28)  = Cj  q + ^49 
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B (1) 


B (2) 


B (3) 


B (4) 


C11D131  + C12°13  + C13°13  + C14D11 


C15  (D134  + D136)  + C16  (D11  + D11  ) 


C17  (D118  + Dlll0)  " C18  (D139  + O 


C11D23  + C12°23  + C13°23  + C14°21 


C15  (D23  + °23  ? + CI6  (D21  + °21  ) 


ft  1 n 911 

C17  (D21  + °21  * " C18  (D23  + °23  * 


C11D331  + C12°33  + C13D33  + C14D31 


C15  (D33  + °33  ) + C16  (D31  + °31  ) 


10 


11, 


I-  C-, (D-jn  + D-,i  ) - C^g  (D33  + Dgg  ) 


'17  '31  31 


= --  2^5^  ^2  cos2  + (Cg  + 
- C3g)  sin2  4>-J  sin  01  cos 


f 2^4^  T32  T33  (Cg  c2  + c37  “ C3g) 
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+ (C34T131  + C35T12  } D23 


" (C34T23  + C35T22  * °13 


+ (C34T13  ” C35T12  * °23 


(C34T23  " C35T22  * D13 


•ft 


B (5)  = ij^2  ^C2  cos2  ^ + (Cg  + C 


37 


- C3g)  sin2  ('-J  sin  0^^  cos  0^ 


- cos  01  [(2  C2  + C38)  cos2  4>1  + (2  Cg  + 2 C 


37 


1 


- C38)  sin  o. 


J 


+ 20^4)^  sin  <j>3  cos  <j>3  (C2  ~ ^6  ~ ^37  ^"38^ 


I 


[C19  (D132t11  + °2 3 T21  + °33  T3i  ) 


+ C20  (D133'ril1  + °2  3 T21  + °3 3 T31  ) 


3m  1 . ^ 3m  1, 


3m  1 


m ^ 1 r\  m > 


+ C21  (D11  T11  + °21  T21  + °31  T31  ) 


4 1 4 1 4 1 

” C25  (D13  T11  + °2 3 T21  + °33  T31 


+ D '''+  D ^ + D 

°13  T11  23  21  33  a31  ' 


D-31 


+ C28  <(DU5  + DU7)  ''ll1 


♦ (D215  + D217)  t217  + {D315  + D317>  T311> 


- CJ2  <fD118  + Dlx  } TX1 


+ {D218  + D2110)  ^l1  + fD^8  + “n1’  ^l1’ 


31  31  J 31 


+ ({D109  + D,,11}  TnX  + (D239  + D23n}  T21] 


'34  ' 13  13 


+ (D339  + °331}  T 3 1 1 >]  COS  *1 

- [C33  ({DU8  " DU10)  Tll1  + {D218  - D2110)  T21l 


+ (d  ® - n 
+ 1D31  u31  * 31  ' 


_ C35  ^°139  " °13  * T11  + *°2  3 “ °23  * T21 

+ ^ - D33U}  T311}]  Sin  ♦! 


B (6)  = (^2  sin2  81  - 6^)  (C2  - Cg  - C37  + D3g)  sin  ^ cos  ^ 
+ 4^0^  cos  8^  j^(2  C2  + C3g)  cos2  <J>^  + (2  Cg  + 2 C37 


- C38)  sin  4> 


J 


_ C19  ^ °13  T12  + °2  3 T22  + °33  T32  ^ 
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B (7)  = (C3  - Cg)  02  sin  02  sin  02  cos  <J>2 


• • 2 *2 

- 2il/202  sin  ^2  cos  6 2 ^C3  COS  ^2  + C5  sin  ^ 2 


2 ’ 5 


+ C39  - C40> 

- 2i2*2  T32  T33  (C3  " C5' 


9 * 2 2 
+ 0 2 4>2  cos  ^2  ^ C3  s*n  4>2  + ^ Cg  cos  (j>2  + C39) 


C (D  ^ - D 

C19  'U13  X2 3 U23  i13  ; 


C22  (D133t23  ‘ D21  T13  ) 


~ C23  (D11  T23  " °21  T13  ] 


+ C26  [(D13  + D13  5 T23  ~ *°23  + °23  * T13  ] 

C27  [(D13  ” °13  * T22  ” (D23  ” °23  * T12  ] 

- C29  j^(D1]L  + D1X  ) T23  - (D21  + D21  ) T13  ] 

+ C30  [(D11  ” D11  ) T22  " {D21  ” °21  ) T12  ] 
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^ 2 

B (8)  - (C,  cos2  *,  9 c5  sin2  *2  9 C39  - C40>  *2  sin  02  cos  02 


'2  5 


- (2  C3  cos2  <I)2  + 2 C5  sin2  <{>2  + C39)  4>2$2  cos  02 
+ 2 (C3  - C5)  624>2  sin  <t>2  cos  <|>2 

+ (C19  (D131t112  + D231t212  4’  D331t312) 


+ C22  <D133t112  + D23JT21*  + D33JT31  1 

+ C23  (D113t112  + D213t213  + ^lSl2’ 

- C26  [(D134  + D136)  Tll2  + <D2  34  + D236>  T212 
+ (D334  * D336)  T31  ] 

+ C29  [(D115  + Dll7’  TH2  + <D2l'  + °217)  1212 


3_  2 . „ 3m  2, 


+ (D315  9 D317)  T312]  }cos  *2 
+ (C,7  [(D134  - d136)  Tu2  9 (d234  - D236)  t 


<D33  " °33  ’ T31  ] 


• C30  [‘“li3  ‘ Dll7)  Tll2  + (D215  ' °217)  T212 


+ (D315  ' D317)  T3i2]  1 Si"  *2 
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B (10)  = 


+ 


+ 


+ 


B (11) 


- (C13p3  + C41  - C42)  [e32T313  sin  *3  cos  i. 


'42 


;t333  (§3 

sin  6 ^ cos  < 

C13p3  + 

2 C41  - C42> 

2 1 

• • 

cos  <{>2 

cos  © 2 ® 3 4>  3 ■ 

_ 1 

3 13, 

(D13  T2 3 

- n T 

U23  A13 

. 1 

3 _ 1 3, 

°13  T21 

- D t 

23  A11  ‘ 

2„ 

3 ^ 2m  3 

(D1  3 T23 

u2 3 A13 

2m 

3 2m  3, 

°13  T21 

_ n t 
U23  A11 

, 3 

3 ^ 3m  3 

°13  T21 

°23  T11 

, 3 

3 3 3 

(D11  T23 

- D T 

U21  A13 

,3P3  + C41  " C42)  COS 

C13p3  + 

2 C - C ) 

41  42' 

. 2 ' 

• • 

sin  4>  ^ 

4>2<!>3  cos  8 ^ 

2 (C13P3  + C41  “ C42)  °3<?>3  sin  (j)^  cos  (j)^ 
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+ [C20  {D131t113  + D231t21  + °33  T31  } 

+ C22  (D132t113  + °23  T21  + D33  T31  } 

+ C24  (D113t113  + D213t213  + D313t313)]  COS  *3 

- C21  [ (D]l31  cos  1^3  + D^1  sin  sin  e3  + D331  cos  03] 

- C23  [(D132  cos  i|»3  + D232  sin  j|»3)  sin  93  + D^2  cos  e3] 

- C24  [(D133  cos  «3  + D233  sin  ^3)  sin  03  + D^3  cos  03] 


2 2 ® 2 

B (12)  = (C13P3  + C41  - C42)  (if 3 cos  93  - 93  ) sin  4,3  cos  ^3 


C13p 3 + 2 C41  " C42J 

cos  ({>3 

2 1 • • 

sin  <f>3  ^36 

3 COS  9 3 

1 3 

„ lm  3 

L _ 1_  3. 

(D13  T12  + 

D23  T22 

+ d33  t32 

2 3 

2 3 

„ 2_,  3 

(D13  T12  + 

°2  3 T22 

+ D33  t32 

,dh3ti23  + 

°21  T22 

3 + D313t32 
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B (13)  - [c?  + C43  - C44  + (Cg  - C45 


2 i*2 

+ C4g)  cosz  a5J  04  sin  e4  sin  (f>4  cos  <p4 

* 2 2 2 2 2 
+ Cg04  (sin  04  - cos^04)  sin  i|>4  cos  \j>4  s.'.n  <p^  cos  a,. 

- (Cg  - C45  + C46)  (e42  - j>42)  cos  e4  sin  <p4  sin  «5  cos  a5 

- { [c7  + C43  - C44  + (Cg  - C45)  cos  a5 

- C46  sin2  a5]  cos2  <p4  - (Cg  - C45)  sin2  a5 

" c46  sin2  <P4  + c46  cos2  a5  > 2 ^464  sin  04  cos  04 
+ 2 (Cg  - c45  + c46)  ii/404  (COS2  e4 
2 

- sin  ©4 ) cos  <J>4  sin  a5  cos  a5 

- 2 [c?  + C43  - C44  + (Cg  - C45 

+ C46)  cos2  ct5]  4 4> 4 cos2  04  sin  <j>4  cos  <j>4 

- 2 (Cg  - C45  + Cec.)y4l4  sin  e4  cos  e4  sin  <f,4  sin  a4  cos  a5 

r 2 2 *1  2 2 

+ 2 { [c8  sin2  e4  + (C45  - C46)  cos  04J  cos  <j>4  - Cg  cos  q4 

+ Cg  sin  <f4  - (C43  - C4g)  sin  04  ) 1^40:5  sin  05  cos  05 

+ 2 (C8  - C45  + C46)  i‘4 a5  (cos2a5 
2 

- sin  ctg)  sm  04  cos  04  cos  <j>4 


D-39 


+ s*n  ^4  COS  ^ 4 s*n  ®4  cos  ®4  s^n  $4  008  $4  cos^<^e 

+ 644>4  cos  e4  {(2  C?  + 2 C43  - C44)  sin2<|)4  + C44  cos2^ 


+ D2  C8  ” C45  + 2 C46^  Sin  ^4  + C45 


2 

cos  a. 


f 


m 


w 


n 


+ C45  sin2a5> 

• * , # 2 
” -•  cg84a5  sin  ^4  cos  ^4  s^-n  64  cos  64  sin  <f> 4 sin  CX5  cos 


- 04a5  sin  e4  sin  <pA  (2  CQ  - CAQ)  cos  ac  + (C 


8 45 


'45 


K 


2 1 

” 2 ^4g)  s^n  a5j 


+ 2 (C8  “ C45  + C46)  ®4^5  cos  64  sin  ^4  cos  $4  sin  a5  cos  <*5 
+ ^4«5  cos  0 4 cos  <j>4  [(2  Cg  - C45)  sin2  a5  + (C45 


- 2 c46>  cos 


“5! 


" 2 (C8  " C45  + ^4^5  sin  e4  sin  a5  cos  a5 
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9 50 


*8T33  > *8 


*“  C34  (D13  T11  + °23  T21  + D331<r319) 


*"  C35  (D12  Tll9  + °22  T21  + D321t319) 


C (D  ^ + D ^ + D 

°36  'Ull  A11  + U21  i21  + U31  i31  ' 


(C10  ~ C49  " C50}  ®102  Sin  01O  3in  *10  C°S  *10  sin2aH 
<C10  + C49  • C50>  (8102  - +102>  T3210  sin  “ll  cos  all 


h 2 | C9  ~ C47  + C48  " [<C10  ” C5o’  cos  *10  + C49.  sln  a 


11 


h (C10  + C49  COS  *10  " c50)  COS  a 


11 


^ C49  sin2*10|  J.1O01O  sin  910  cos  e10 


h 2 (cio  - C49  " C50)  *10910  (GOS  6 


10 


2 

sin  9^q ) cos  4»^q  sin  cos 


2 (C10  + C49  " C50)  01O^1O  T33  T32  Sin  all 
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+ 2 j [C10  + (C49  ' C50>  cos2*lo]  °Os20lO 


- (C,0  oos24.10  + C49  - C5Q)  sin20lo 


- C sin  410  *10nn  sln  “11  cos  “11 


+ 2 (C1Q  + C49  - C50)  *10in  l=os  an 
o 10 

- sin  au)  T33  sin  910 

+ jc47  + ['2  cio  + 2 C49  C50)  Sin  h0 

2 1-2 

+ c50  cos  hoj  Sln  all 


+ C50  cos  au[  01O*1O  cos  el0 

+ | [(2C10  + 2 C47  - C50)  sin2ail  + C5Q  cos2ail]  sin  e10 
- 2 <Cl0  + C49  - C5Q)  T3310  sin  au  cos  ccjJ  «lo“U  *ln*i0 


2 (c10  + C49  ' C50l  sin  610  Sin  “li  C°S  “11 


- [(2  cio  + 2 C49  ‘ C50>  cos2“ll 


+ C50  sin 


in2au]  T 


10  . . 

33  hoall 


♦ C32  ‘“uSl10  - DH  V°>  - c33  <D121t2110 


D221tU10) 
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B (27)  = (C1Q  + C49  - C5Q)  ^102T3210  sin  e1Q  sin  «n  cos  au 


+ (C10  + C49  - C5Q)  (^10  cos  0 


10 


“ ) sin  cos  <J>^q  sin 


+ 2 (C1Q  + C4g  - C5Q)  (01O  cos  <j> 


10 


+ a13)  i|»^0  sin  01Q  sin  cos 


+ C4^  + £(2  C^Q  + 2 C49  - C^q)  cos  (j> 


10 


4 C50  sin\oj  sin2“n  4 c50  00s2c,nl  *ioeio  cos  eio 


+ [<2  C10  4 2 C49  C50*  sln  “ll  4 C50  COS  “ll]  (*10T33 


®io  sin  ^lO*  all 


• • 


2 (C10  + C49  " C50)  ^lO®!!  Sin  all  COS  all 


- [C34  (D131t1210  + D231t2210  4 D331t3210> 


C (D  + D ^ + D 

^35  u12  a12  u22  1 22  u32  A32  ' 


- C36  (D1110t1210  4 D2110t2210 


0 10m  10.1 

°31  T32  U Sln  all 
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B (28) 


(C10  + C49  ‘ C50}  ho  ^sin2aH  " "^“ll*  T3310  sin  610 


(C10  "Am  + C49  - C50?  Bi*\o 


'10  49  50' 


[c10  + <C49  - C50;  cosJ»lo] 


2„ 

cos  0 


10 


I-  C10  sin2((.10  i^102  Sin  au  cos  an 


f (C10  + C49  ■ C5o'J  " + 10  sln  +10 


10 


+ 2 + i o (0ior32‘  oos  no  no 


+io  • +io  sin  °io)]  sin  an  cos  an 


* [(2C10  - 2C49  - C50>  00s2all 


+ C5Q  sin2au 


] +10+10  s-‘;  910  Sin  +10 


<2C10  + 2C49  - C50*  sin2“ll  + C50  cos2“ll]  <®.’«  sin  +' 


'10  T10 


” *10T3310)  +10 

♦ = 34  '“nV1  + D231t21U  + “asSl^ 

- =35  <D12lTllU  + D221t21U  + ^Sl^ 


- =36  <DJl10Tlin  + D2110t21U  + D3110t31U> 
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Below  are  defined  the  constants  D^n  used  in  the  B (N)  on 
the  previous  pages, 
n = 3,  8,  10 

Dnn  - <*„2  + 5n2)  Tun  - 2 *n9n  sin  *„  sin  6n 


D2in  = (*n2  + O T2l"  + 2 cos  *„  sin  9n 


°3in  - - 9n2  sin  9n 
n = 1,2 

n O O Y\  1 

= (to  + to  j T.  0 +0  cos  to  sin  0 sin  <b 

12  'Yn  yn  12  n rn  n Yn 


+ 2 10^  sin  to  cos  0„  sin  <J> 
n n n n n 


+ 2 Vn  T23 


n 


- 2 0 to  cos  to  cos  0 cos  <t> 
n n n n n 


D„„n  = (4 » 2 + 4>  2)  T00n  + 0 2 sin  to  sin  0 sin  to 
22  rn  ^n  22  n rn  n n 


2 to  0 cos  to  cos  0 sin  to 
n n rn  n rn 


- 2 | T.  _ 

rnYn  13 


n 


- 2 0 to  sin  to  cos  0 cos  to 
n n n n n 


D = (0  2 + | 2)  T„n  + 2 0 to  sin  0 cos  <{> 
32  n Tn  32  nyn  n yn 
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n = 1/  2,  3,  4,  6 


°13n  = <*nZ  + *n]  T13n  + *n  COS  K sin  6n  OOS  *n 


+ 2 ^n6n  sin  4n  cos  0n  cos  4n 


" 2 Vn  T22 


+ 204  cos  4 cos  0 sin  4> 
nTn  rn  n n 


D23n  = (i,n2  + O + C sin  K sin  0n  cos  *. 


-240  cos  4 cos  0 cos  <b 
yn  n rn  n Tn 


+ 2 Vn  T12 


+204  sin  4 cos  0„  sin  4„ 
nTn  rn  n n 


D33n  * (5n2  + O ' 2 Vn  sin  6n  Sin  *n 


l = 5,7;  n =•  l - 1 


DH*  = + ^2)  Tll£ 


+ 0 cos  4 ( cos  0 sin  a„  - sin  0n  cos  4 cos  a0) 

n n n X/  n n x. 


• 2m  n 

' 9n  T11  COS  a£ 


+ 2 4 0 T0,  ^ sin  4 

rn  n 31  n 
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+ 2 Vi,  T22n  co^ 


+ 2 Vt  T23 


n 


2 Vn  T32  cos  ♦n  cos  “t 


2 Vl  T33  005  *n 


n 


+ 2 a T10  sin  a 
n 12 


r,  l , • 2 , • 2.  m i 

D21  + ai  ) T21 


+ 0 sin  ip  (cos  6 sin  a„  - sin  0 cos  <p  cos  a.) 
n n n H n n X. 


' *n2  T23n  cos  aH 


- 2 iji  0 T_.  cos  ij; 

n n 31  rn 


2 Vn  T12n  cos 


2 V*  T13 


n 


2 Vn  T32  sin  *n  cos  at 


- 2 0 a„  T__  sin  ^ 
n i 33  Yn 


+ 2 Vn  T22n  sln 


D-G9 


D3l‘  * <«n2  + “/>  T3l" 


' ♦n2  T33n  cos  aH 


+ 204  sin  0 sin  <f>  cos  a0 
n n n n x 


+ 2 0 a.  (cos  0 cos  a0  + sin  0 cos  4 sin  a.) 
n x n Jt  n n x 


n 


+ 2 Vl  T32  Sin  a 


m=9,ll;n-m-l 


_ m / ; 2 • 2.  m 

°13  “ (iJn  am  ) T13 


- 0 2 cos  4 (cos  0 cos  a - sin  0^  cos  ® sin  a ) 
n n n in  n n m 


• 2 _ n 
+4  T. _ sina 
n 13 


m 


• • „ m 


+ 240  T_ _ sin  ip 
rn  n 33 


m 


n 


-244  T_  ^ sin  a, 
rnyn  22 


m 


+ 2 i a T01 
n m 21 


in 


+ 2 5 (j>  T0On  cos  4 sin  a 
nyn  32  rn  m 


- 2 Vm  T3l”  COS  ♦»  + 2 V«  T12n  cos 


m 
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“as"  * <^n2  + %>  *23" 


- 0 sin  \b  (cos  0 cos  a - sin  0 cos  <J>  sin  a ) 

« "n  n T n tn 


n 


n 


n 


rn  m 


• 2 „ n 
+ d>  sin 

Yn  23  m 


i 


m 


- 2 Vn  T33  cos  ♦n 


■ 


n 


+ 2 H sin  a 

rnTn  12 


m 


% i 

f t 

[ 

l I 

I.  I 

f : 


• • m 

- 2 *n“m  TU 


+ 2 Vn  T32n  sin  *n  sin  am 


m 


- 2 0 a T sin  to 
n m 31  n 


I 


E ' 


n 


+ 2 <b  am  T_0  cos  am 
Yn  m 22  m 


I 


D33m  “ <C  + “rn'’  ^ 


• 2 n 

+ <b  T00  sm  a 
Yn  33  m 


-204  sin  0 sin  <j>  sin  am 
nyn  n rn  m 


+ 2 0 ct  (cos  0_  sin  + sin  0 cos  6 cos  a ) 
n in  n m n n m 


n 


+ 2 <b  a T „ cos  a 
Yn  m 32 


m 
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APPENDIX  Fi 

ELEMENTS  OF  VECTOR  {P} 


E-l 


0 


F =0 
P2 


“Cl* 


F =0 
P4 


F = C, , g sin  0 , cos  rf, , 

P5  11  1 T1 


= C„gT 


Pe  32 


F =0 


F = C,-g  sin  0 _ cos  A _ 
Pg  LZ  Z Z 


F = C,0gT  ' 
P9  12  32 


F =0 
P10 


Fp  = (C13  sin  03  cos  ,j,3  + Cl4  cos  q^)  g 


11 


F = C, ,gT0  ' 
p^2  133  32 


F =0 
p13 


14 


= -g  [c15 


cos 

^4 

cos  a 

cos 

64 

sin  a 

Fp  - - <c15  + C16  008  “5>  ^32 


'15 


Fp16  = ' C“sT« 


F =0 
P17 


Fp  = "g  C15  [sin  06  COS  *6  + C16  (sin  6 6 C0S  ^6  COS  a7 
1 8 


- cos  0g  sin 


Fp  = ~ (C15  + C16  COS  a7>  gT32 


'19 


F = -C  figT 
p20  16  33 


F «:  0 

P21 


Fp  = g cos  08  T C18  ( cos  0g  cos  a,. 


- sin  0g  cos  <j)g  sin  ag)j 


Fp23  = -C189T32  sln  °9 


FP24  Cl8gT3l' 


p =0 

P25 


F = 9 [c17  90S  ei0  + C18  (COS  010  CCS  au 
P26  L -I 

ho  sin  “ll’J 


- sin  e10  cos 


E-3 


Bf 


I 


I. 


1 


1 


i 


I 


I' 

5 


27 


10 

= “C18gT32  Sin  all 


p28 


~ C18gT31 


11 


E-4 


i 


1 


i 

i 


APPENDIX  F 

ELEMENTS  OF  VECTOR  {R} 


F-l 


Joint  Resistance 


Joint  Identification 

1 : spine,  between  segments  1 and  2 (between  T-12  and 

L-l) 

2 : neck  (between  C-7  and  T-l) 

3 : right  shoulder 

4 : right  elbow 

5 : left  shoulder 

6 : left  elbow 

7 : right  hip 

8 : right  knee 

9 : left  hip 
10  left  knee 

Elements  of  {R}: 


10  . 38 

r3  = ■ i?  1 (Mi + w 2 


28) 


where  81  cos  (T13  T13  + T23  T23  + T33  T33  ) 

= cos  ^ H, 


F-2 


for  i=  1,  2,  3,  5,  7,  and  9 are  defined  below 
H1  = T13lT132  + T23lT232  + hiS* 


H2  “ T13  T13  + T23  T23  + T33  T33 


H3  T13  T13  + T23  T2 3 + T33  T33 


H =T  ^+T  6+T  ^ 

3 13  x13  + A23  L22  133  L 22 


H = T ^ + T ^ + T ® 

7 11  X11  + l21  21  31  i31 


H = T XT  + T + T 

H9  T11  T11  T21  T21  T31  T31 


whereas,  for  the  elbow  and  knee  joints, 

B4  = “5 
S6  = a7 
B8  = a9 
B10  = aU 


The  derivatives  required  for  computation  of  are 


1 3H. 


aqj  Vi  - hi2  9qj 


F-3 


and 


Si 


H. 


i1  - hi 


(i  — 1 / 2 , 

3 , 5,  7, 

9, 

• 

Defining 

for  i = 

1, 

• 1 2 

• 12 

61-T13T13  + 

T23  T23 

+ ' 

1*  2 

m 1*  2 

+T13T13  + 

T23  T2 3 

+ 

• 2 3 

• 2m  3 

G2  = T13  T13  + 

T23  T23 

+ 

2*  3 

m 2*  3 

+T13T13  + 

T23  T23 

+ 

• 2 4 

G3  = T13  T13  + 

T 2T  ^ 
i2  3 23 

+ 

2*  4 

m 3*  4 

X13  13 

T23  T23 

+ 

• 2 6 

• 2_  6 

S = T13  T13  + 

T23  T23 

+ 

2*  6 

m 2 * 6 

+ i13  13 

T23  T23 

+ ri 

• 1 8 

•18 

G7  = T11TH  + 

T21  T21 

+ 

1*  8 

m 1*  8 

+ T T + 

+ A11  11 

T21  T2 1 

+ 

• 1 10 

• 1_ 

10 

G9  = T11  T11 

+ T21  X21 

1*  10 

_ 1- 

10 

+ T11  T11 

+ ^21  l21 

1„  2 


, H 2 

33  33 

■ 2m  3 


. 2t  3 
33  33 

; 2m  4 


, 2 • 4 

1 T 
33  33 

, 2„  6 


2*  6 
1 Aip 

33  33 

•18 


, 1-  8 
T 

31  1 31 


, iT  10 

31  31 

, 4 10 

•31  *31 


F-4 


where  the  time  derivatives  of  the  transformation  matrix 
elements  are  given  by 

*Un  - - *nT2in  + ®nT3in  008  *n 


• n • n • n 

T„,  = ip  T " + 0 T “ sin  & 

21  Yn  11  n 31  n 


T_, n * - 9 cos  0 
31  n n 


*13"  ■ ' *nT23n  + 008  *n  ' 


*23n  ' Vu  + 5„T2in  808  ♦»  ' *nT22n 


*33"  - SnT3l"  808  *n  ' *nT32° 


The  non-zero  elements  of  {R^  are 


“ r=  (“l  ' 7=4  (T131t232  - T231t132) 

+ r==/"l  - 4=)  (T111t218  - T211tv.8) 

+ 7=H(M9 ' J^=%)  (TlllT2ll° ' T2i1tu1°) 

V-l  - H 2 V V1  - V 


P-5 


R5  = 


7===.  fMi  - 1====)  <T111t132  + T211,] 

11  13  21 


23 


+ 'I311t332)  COS  h 


r 

Vl  “ H?2  V Vl  - k72 


(Til  COS  ^1 


8 8 1 
+ T21  sin  ^i>  sin  6i  + T?i  cos  9i 


(Tn10  cos  ^1 


R6  " 


+ T2i^  sin  ^i>  sin  9i  + T3iA°  cos  ®i] 

= (Ml  ’vrrtj)  ^n^132  + 122X1232  + T32lT33^ 


r? = (Mi  * ^is2t231  ■ 1232113 ' 
+ ■■  ■ --  1 _ (m2  - - 2 - ~ \ (t 

\ y 1 - h22/  ' 


2 3 2 3' 

13  x2 3 a23  1 13 


+ L__  /„  - J5  G5  \ /T  2 6 . T 2 6\ 

I y lM5  1 j)  ri3  X23  L23  i13  ] 

yi  - v ^1  - h5v 


F-6 


ML  *** 


+ T313t332)  COS  *3 


R12  = ‘ f==  (M2  " 7==)  (T123t132  + T223t232  + T32  r33  ) 

- h22  v yi  - V7  x 

,i!  'rt5  ("■  ■ ',,v) 

r14  - (»3  - -p=)  (T114t132  + *2lV 

yi  - h32  v tji  - h3  / \ 

+ T334t332)  °OS  *4 

*15  = - - 7==)  (T124t132  + T224t232  + T324t332) 

V1  - H,2  V ^ - h32/  ' 


R16  “ " (M4  + J4  a5} 


1 /„  J5  °5  \ fT  6 2 6 2 

= 7 =7  ( 5 _ r ===)  lTll  T13  + T21  T23 

Vi  - h52  Vi  - «52/  X 

+ T316t332)  cos  "“6 


+ T^1  sin  ^Q)  sin  0g  + T^1  cos  9 


+ T^1  sin  i|»10)  sin  e1Q  + T31X  cos  e10 


'10' 


10  31 


'28 


" " ^M10  + Jio  an^ 


F-9 


APPENDIX  G 

ELEMENTS  OF  VECTOR  {Q} 


G-l 


*1  t 

i 


°e  ■ \ [\  'll1  - {yz1  - P1>  ’l/] 

+ \[y*l Tsl1 ' (Yzl ' pl>  Tz2^ 
+ Pzl  [Yyl ^ ' (Yzl  ■ pl>  "32^ 


[■  ,f*2  + 


F*  > *12*  - (Fy2  + Fy3>  *22" 


- (Fz2  + %>  T321] 


+ Pi  T(Fx  + F + F + 
1 L 8 X9  X10 


T ) T - 
X11  12 


+ (F  + F + F + F„  ) T,, 

y$  y9  yio  yn  22 

+ (F  + F + F,  + F,  ) T..1 

z8  z9  z10  Z11  32  J 


f(F  + F 

L x8  X 


“ F - Fx  > T13 
9 x10  11 


+ (F  + F - F - F„  ) Too 

y8  y9  yio  yu  23 

+ (F  + F - F - F ) T,,1] 

z8  z9  z10  Z11  33  J 


G-4 


‘ vm$l ■-' 


<>10  - Px3  [-  <%  + V V - S *22*  - 0,3  + P3>  T233] 

+ Py3  [<Yx3  + V TU3  + Yy3  T123  + <Yz3  + p3»  T133] 

«U  ’ % [<Yx3  + V T313  + Yy3  T323  + *%  + P3>  T333]  » 


COS  4'. 


+ % [<y*: 


F„  r<Y„  + V T313  + yy  T323  + <V  + P3>  'r333|  sln  ♦: 


z3  3 33 


- F Ry^  + Lc)  cos  e3  + Yy^  sin  63  sin  <l>3 
+ (Yz  + P3)  sin  63  cos  <I>3J 

<>12  ■ Px3  [V,  T133  - <%  + P3>  T123] 

+ Py3  [Yy3  T233  ' (\  + P3>  T223] 

+ % [Yy3  T333  ‘ ,Y.j  + P3>  T323] 

°21  ■ Px8  [-  <Yx,  + P8>  *2l'  - Yy8  T228  - Yz„  T238] 

Py8  [«\  + V Tll8  + Yye  T128  + Yz8  T138] 

♦ Px,  [-  L8  T218  - Yx,  *21*  - \ T228  * <%  ' V T239] 

+ V R "ll®  + \ TU9  + Yy9  T128  + ,Yz9  - V t139] 


i 


G-6 


Q22  Fx, 


[<Yx8  + P8>  T318  + Yy„  T328  + \ T338]  c0S  *8 

+ Fy8  [(Yx8  + Ps)  T3J-8  + Yys 1328  + Y*s T33^  sin  4,8 


- F 


[(YX8  + p8>  cos  8 8 + Yy8  sin  8 8 sin  *8 


+ yz8  sin  8 8 005  *8j 


+ % [s  T318  + Yx,  T319  + Yy9  T328 


+ <Yx9  - V T33  J 003  *8 


+ F 


[L8  T31  + Yx0  T319  + YyQ  T32 


+ <\  " L9)  T33  ] sin  ^8 


- Fz  ^Lg  cos  9g  + yx  (cos  0g  sin  otg  + sin0  g cos  <J>{ 


+ y„  sin  0g  sin  <j>8  + (yz  - Lg)  (sin  0H  cos  <j)Q  sin  c 


8 T8 


- cos  0.g  cos 


“s’] 


°23  * % <Yy8  T138  - Yz8  T!28» 


+ F (Y  T 8 - y T 8) 

y8  y8  23  Zg  22 


cos  ag) 


9 


G-7 


,1 


+ F.  T 8 - Yz  T 28) 

Zg  Yg  33  Zg  51 


+ % ['  % Tl28  C0S  0,9  + Yy9  TnS  ‘ N ’ ^ Tl28  Sl" 

+ Fy9  ['Y*9  T22*  cos  “9  + Yy9  1238  ’ ‘N  " L9>  T22'  sin  ^ 

+ Fz(J  [-%  T328  003  a9  + Yy9  T338  ' (YZg  ‘ V T328  sin  “9] 

°24  ' % [-  % T138  + - V Tll9] 

+ Fy,  ['Yx9  T339  + (Yz,  ' V T319] 

+ % ['%  T”9  + <Yz9  ' V 1319 ^ 

Q25  ■ - Fxw  [(Yx10  + PB>  T2110  + Yy10  T2210  + Yz10 


'24  x( 


+ F f(Y  t Pa)  Tlt  ° + Y T , + Y T 10] 

y 10  L xio  8 11  yio  12  zio  Ai3  J 


[l8t; 


xu  L“8*21 


+ Y-U  T21U  + Yya  t2210  + <V  - V T«U] 


+ Fyn  h Tu10  + V *u“  * Yyu  ’u” 
+ (Yz,,  - V T13U] 


G-8 


°26  ' Fx 


10 


[<\o  + P8>  T3110  + "y10  T3210  + T3310]  “s  *10 


m 


+ \o  [So  + P8>  ^ + Yyi0  + \o  T33X0]  si"  *10 


m. 


- F 


'10 


[<yx10  + p8>  oos  °10  + Yy,„  sin  610  sin  *10 


10 


&- 

§ 

i 


Wi 


+ Y .0  d>  1 

Z10  sin  10  cos  ioJ 


+ (PXll  cos  *10  + Fyu  sin  *10>  [L8  T3110  + yxn  T31U 


I 


+ yyuT3210+  <yz„  ' V 


'11 


T33U] 


- F TLq  cos  0,.  + y (cos  0,_  sin  a., 
z,,  [_  8 10  x,  n 10  11 


'11 


'll 


+ sin  610  cos  *10  cos  au>  + Yy^  sin  610  sin  *10 


1 

I 


£ 


+ (yz  - Lg)  (sin  01O  cos  <{>10  sin  a ^ - cos  01Q  cos  a^)] 


10  11  J 


°27  = Fx  (Yv  T1310  “ T.  5 10 ) 

X10  y10  13  12 


'1C 


+ F (v  T - v T 

*10  y10  23  Yz10  22  ’ 


+ Fz  <V  T3310  - yz  T3210) 
10  yl0  ' n 


'10 


+ Fx  \yv  T1310  “ (Yx  cos  al] 

X11  L yll  ■LJ  X11  11 


G-9 


, v:  £.;/ ' fv 

A- 


>i 


+ Wzn  ’ 1,91  sin  oii>  Ti2l0i 

+ P |>„  T,,10  - <YV  COS  a-. 

yu  L yn  23  n 

+ (Y  - V sin  “u>  T2210] 


+ F*n  [Yyu  133  ' Ni  005  “u 

+ {yzu  ‘ V Sin  011>  T3210l 

°28  ' V [%  TnU  + <V  - V TnUl 


*JJU  + (yz„  - V *21U] 


yll  I xu 


['V  ■ V 
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H.l  CALCULATION  OF  MOMENT  OF  INERTIA  FOR  SEGMENT  2 


.. L 

t 3.04 


7.18 


0 


Combined 

Shoulder 

Mass 


Thorax 


A 

1 r 

5.15 

Upper  Torso 
(Segment  2) 

10. 

v 

For  Dempster's  segments,  idealized  above  (a)  the  combined 
mass  of  the  shoulder  segments,  mg  = 0.1054  M where  M is 
total  body  mass.  The  mass  of  the  thorax  segment, 
mfc  = 0.1097  M.  For  the  50th  percentile  male  with  total 
body  weight  = 161.5  lb, 


m = 0.1054  (161.5/386.4)  = 0.0441  lb-secI 2/in 
mt  - 0.1097  (161.5/386.4)  = 0.0459  lb-sec2/in 

Dempster's  moments  of  inertia  for  transverse  axes  passing 
through  the  mass  centers,  adjusted  for  weight  of  161.5 
lb  are 

I = 2 (0.5613)  = 1 .1226  lb-sec2-in 
s 

= 1.1149  lb-sec2-in 


Using  the  parallel-axis  theorem, 

Iy2  = 1.1226  + 0.0441  (5.15  - 3 . 04 ) 2 
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H.2  CALCULATION  OF  I AND  I 

xn  2n 

1.  Approximate  segments  1,  2,  and  3 by  ellipsoids 


x 


Segment  1 


T - m /H2  x „2\ 

Xx  " 5 (b  + c } 


t = § (a2  + c2) 


t m , 2 , ,2. 
1z  ■ 5 (a  + b 1 


2 

I = 4.331  lb-sec  -in 
yl 

Assume 

a = c = (buttock  radius) 

b = L_.  + R..  (1/2  hip 
11  breadth) 

For  50th  percentile 

a = c = 4.0f  b = 13.9/2 


2 

Xx  = f [(^r)  + (4)2]  = 64.3  m/5  = I, 

I = § [(4) 2 + (4) 2]  = 32.0  m/5 


2 

I = I = 2.009  I = 8.703  lb-sec  -in 

Xx  zi  yi 
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Segment  3 


T 


y3 

a 


0.311  lb-sec2-in 


R^  = 4.0  in  = b 


c = 


p^  = 5.8  in 


l 

I 

Ix  = ? [(4)2  + <5-8)2]  = 49.6  m/5 
3 

I = I 
y3  X3 

Iz  = f [(4) 2 + (4) 2]  = 32  m/5 

J 

I = 0.311  lb-sec2-in 
x3 

I_  = 0.6446  I = 0.2006  lb-sec2-in 
z3  y3 
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H-7 


vim^nx S^rc 


m i 


Forearm  and  Hand  Segments 


I.  = 0.218  lb-sec  -in  = I 


Y = R5  1.85  ; l = L5  = 13.23 


I = I =1/4  mR I  2 + 1/12  mLc2 

y5  z5  5 5 


I = 1/2  mR  2 = 0.0241  lb-sec2-in 

x5  •* 


I = I 
x?  x5 


S = % " ry5  ''  \ 


Thighs 


I = 1.270  lb-sec  -in  = I 

y«  zi 


Y = Rg  = 3.55  ; £ = Lg  = 16.58 


I = 1/4  mR  2 + 1/12  mL  2 
z8  8 8 


I = 1/2  mR 


8 


8 


= 0.307  lb-sec2-in 


Xx  " 

x10  x8 


I = I =1=1 


10 


J10  7 8 


ya  zs 
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'•  ''"•?  ?'-  ’"V  ■>  Wv'*  , - ’-' 


Leg  and  Foot  Segments 


2 . 


I = 1,192  lb-sec  -in 
y9 


y = R9  = 2.3  ; % = Lg  = 17.31 


Iy  - 1/4  mR,2  + 1/12  mLg2  = Ix 


1 = 1/2  mR9‘ 


2 

0.1199  lb-sec  -in 


Ixll  yll  X9  y9 


h = x2 

Z11  z9 
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